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Abstract 

A four part series of lectures on the connection of statistical mechanics and quantum 
field theory. The general principles relating statistical mechanics and the path integral 
formulation of quantum field theory are presented in the first lecture. These principles are 
then illustrated in lecture 2 by a presentation of the theory of the Ising model for H = 0, 
where both the homogeneous and randomly inhomogeneous models are treated and the 
scaling theory and the relation with Fredholm determinants and Painleve equations is 
presented. In lecture 3 we consider the Ising model with H ^ 0, where the relation with 
gauge theory is used to discuss the phenomenon of confinement. We conclude in the 
last lecture with a discussion of quantum spin diffusion in one dimensional chains and a 
presentation of the chiral Potts model which illustrates the physical effects that can occur 
when the Euclidean and Minkowski regions are not connected by an analytic continuation. 
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Part I. Introduction 



Abstract 

The connection between the Euclidean path integral formulation of quantum field 
theory and classical statistical mechanics is surveyed in terms of the theory of critical 
phenomena and the concept of renormalization. Quantum statistical mechanics is surveyed 
with an emphasis on diffusive phenomena. The particle interpretation of quantum field 
theory is discussed in the context of nonperturbative statistical mechanics. 

1. Introduction 

There is probably a theorem somewhere in The Critique of Pure Reason which says 
that it is impossible to use language to communicate an idea from one person to the other 
unless the second person already knows what the first person is talking about. Conse- 
quently it is probably logically impossible to teach anyone anything. 

The practical consequences of this unhappy theorem are all too evident. If I wish to 
be clear I must go to great lengths to define and explain my terms. But by the time I have 
achieved utmost clarity I have lost my audience and there is no one left to hear what new 
things I have to say. If on the other hand I wish to keep my audience, I must assume that 
they will use words the same way that I do and proceed directly to my results. But the 
unhappy fact is that since no two people use words to mean the same thing my conclusion 
will be lost on the audience which spends its time in trying to fit my conclusions into their 
definitions of the words used. 

This basic difficulty in the theory of knowledge is exceptionally evident in the presen- 
tation of the relationship between quantum field theory and statistical mechanics. This is 
because the two subjects have developed from very different roots and have very different 
languages and assumptions. 

Historically quantum field theory (QFT) is intimately linked with the classical field 
theory of electromagnetism and with particle physics. Experimentally QFT has been 
intimately connected to high energy physics experiments at accelerators starting with 
cyclotrons in the 30's and 40's and proceeding up to SLAC, CERN, Fermi Lab, KEK, 
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DESY and CEBAF. Theoretically QFT has gone from QED to QCD, the standard model 
and string theory. 

The origins of statistical mechanics (SM) are totally and completely different. Histor- 
ically SM is linked with the theory of heat, irreversibility, and the kinetic theory of gases. 
Experimentally SM has been intimately connected with calorimeters, specific heats, mag- 
netic order parameters, phase transitions, and diffusion. Theoretically SM has gone from 
ideal gases, series (virial) expansions, and Landau theory to solvable models, Yang-Baxter 
equations and holonomic sets of equations for correlation functions. 

It would thus seem from the above that QFT and SM have nothing to do with each 
other. Indeed in the practical functional sense all physics departments do indeed act in 
this fashion. When categories for recruitment are set up a sharp distinction between a field 
theoretician and a statistical mechanician is always made and it is a career determining 
decision on the part of a young person whether QFT or SM is chosen. 

Against this background of almost total separation I have been given the job of ex- 
plaining in these lectures why the two subjects are related. It is not exactly clear what 
the organizers had in mind when assigning me this task. Perhaps the first answer that 
comes to mind is that the path integral formulation of Euclidean QFT is concerned with 
the computation of path integral averages that formally look like 



< 0j >= wy wt , Ze = [ me s^ (L1) 

Ze J 

where f[d<j)] stands for a functional integral over a set of Euclidean quantum fields in d 
spatial dimensions where the distance s is s 2 = x\ + . . . + x% Oj is constructed from those 
fields, and Se{4>) is a function of the fields known as the action which is expressed in terms 
of a local Lagrangian 

S E {(j>) = J dx 1 ...dx d L{ ( j)(x)). (1.2) 

In classical statistical mechanics thermal expectation values are computed as 

V n. P -E(<t>)/kT 
< Oj >= ^ states 3 , Z = e~ E ^/ kT (1.3) 



all states 



where Eaii states i s a sum over au states of the system described by the random variables </>, 
Oj is constructed from 0, E(<j>) is a classical energy functional of 4> , T is the temperature 
and k is Boltzmann's constant. These equations are "formally the same" and thus it may 
be expected that there is a relation between the subjects. 
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Of course and (|1.3| ) are not the only starting points possible. I could, for 

example, choose to start with a Minkowski QFT in d — 1 spatial and one time dimension 
defined by 

< o, >= !MQf™\ Zm = / W]e «.c«/» (L4) 

Zm j 
where the metric is s 2 = x\ + . . . + x 2 d _ 1 — t 2 

Sm{4>) = j dt J dx! . . .dxd-iL(4>(x,t)) (1.5) 
and the quantum thermal average 

< 0j >= , Z = Tre~ H ' kT (1.6) 

where H is a quantum Hamiltonian. Thus from the outset I will recognize that there is 
some flexibility in meaning of the words in the title of these lectures. 



But if it is the relation of (|1 . 1|) to ( |1.3| ) which the organizers had in mind then I am 
in somewhat of a peculiar position because one of the lecturers at this school is Claude 
Itzykson whose recent book Statistical Field Theory [|1J is an outstanding discussion of this 
very relation in a much more detailed and extensive format than I can possibly match in 
these four lectures. Indeed, in the last 10 years many books have appeared [[i | -fTj which 



exploit the relation between ( |1 . 1| ) and ( |1.3| ) . Their titles range from gauge theory to string 



theory to critical phenomena and each has a distinct point of view. The sheer number of 
these books demonstrates the unity of the subjects, a unity which is emphasized by the 
fact that almost all these books, regardless of title or orientation of the author, study the 
two dimensional Ising model as one of the first examples of the subject. 

These four lectures cannot hope to be a substitute for the one semester or one year 
course which can be taught from the existing literature. Therefore I must be selective in 
my choice of topics. 

In this first lecture I will discuss the general principles which lead to the unity of 
( |1.1|) and (|1.3[) . I will approach this from the point of view of the skeptic who needs to be 



convinced. This will hopefully overcome some of the language problems that sometimes 
obscure understanding. 

In the second and third lectures I will illustrate these general principles by presenting 
the field theory construction of simplest of all systems: the Ising model in two dimensions. 
This subject is extremely rich and will illustrate in terms of exact computations practically 
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all of the general phenomena. In lecture two I will discuss the Ising model in zero magnetic 
field and summarize the field theory insights obtained these past 50 years. 

In lecture three I will consider the Ising model in a magnetic field. In this case I will 
present an exact mapping to the simplest model of a two dimensional gauge field coupled 
to a self interacting scalar field. I will then summarize the analytic computations done on 
this model in the last 20 years. We will see that even for this most simple of gauge theories 
that there are many interesting unresolved problems. 

Finally in lecture four the considerations are extended to quantum statistical mechan- 
ics where we deal directly with Minkowski space. We discuss quantum spin diffusion and 
conclude with a discussion of the chiral Potts model as an example of a system where level 
crossing causes the physics in the Euclidean and Minkowski regions to be very different. 



2. Generalities 

In order to be clear it is perhaps desirable to make explicit some tacit assumptions 
that have been made in writing down ( |1 . 1| ) and ( |1.3| ) . 

Perhaps the most basic relation between SM and QFT is not ( |1 . 1| ) and (|1.3| ) but 



rather the fact that both deal with systems in an infinite volume and hence with an infinite 
number of degrees of freedom. A major consequence of this is that the formal definitions 
( |LT|) and (|1.3| ) by themselves have no meaning at all because they are, at best, evaluated 
as — . There is always a further definition needed to make sense of these expressions. In 
the case of SM that definition is embodied in the thermodynamic limit which first evaluates 
( |1.3|) in a finite volume and then takes the limit as the size of the box goes to infinity. In 
the case of QFT the expression ( |1 . 1| ) needs the additional definition which is provided by 
a "renormalization scheme" that usually involves a short distance cutoff as well as a finite 
box. 

In SM the (infrared) thermodynamic limit is treated very explicitly. In QFT the short 
distance (ultraviolet) cutoff is discussed extensively. The difference in focus on infrared 
cutoffs versus ultraviolet cutoffs is often one of the major barriers to communication be- 
tween the two fields and to me seems to constitute a major reason that QFT and SM are 
traditionally considered to be different subjects. 

A second tacit assumption is that ( |TTT| ) and ( |1.3| ) will in the end have an interpretation 
in terms of "particle" or possibly "quasi-particles." The "particle" interpretation of QFT 
is fundamental for its relation to particle or high energy physics and the quasi-particle 
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interpretation of ([O]) is fundamental for such things as the Landau theory of the Fermi 
liquid. This common particle interpretation has the opposite effect from the limiting 
procedures and makes there seem to be great unity between the subjects. Indeed from at 
least the late 50's there has been a great industry in the "application of QFT methods to 
SM" which is ultimately based on this common interpretation. 

There is one further point which is related to both the particle description and the 
infinite degrees of freedom which should be made explicit; namely, the distinction between 
"bare" on the one hand and "renormalized," "collective," or "dressed" on the other. The 
expressions for the fields (p, actions and interaction energies which are used in (|1 . 1|) and 
( |1.3|) before the various limits are taken are often referred to as bare quantities. These are 
the fundamental inputs to the theory and define the problem we are interested in solving. 
It is very common to give particle names to these bare fields such as electron, neutrino, 
quark, or gluon. On the other hand what we observe in measurements are the expectations 
defined by and ( |1.3| ) after the limits have been taken. These expectation values 

are also given names which are often the same at the "bare" fields which appear in the 
expressions for Oj. To distinguish these names apart the words "renormalized," "dressed" 
or "collective" are often used when referring to the expectation values. It is most important 
to realize that because we are dealing with systems with an infinite number of degrees of 
freedom that "bare" and "renormalized" quantities may in fact have nothing in common 
except the name. The bridging of the gap between these two constructs is fundamental to 
the understanding of the connection between QFT and SM. To attempt to make a proof 
for a "renormalized" quantity by first making a proof for the "bare" quantity and then 
adding the word "renormalized" to it runs the grave danger of reducing theoretical physics 
to process of "proof by means of making a bad pun." 

By now the difficulty in communication has become abundantly clear. According to 
your taste I have either 1) taken a large number of well known concepts and tried to make 
them mysterious by putting quotes around them or 2) I have introduced a large number 
of totally undefined words and pretended that they have an unambiguous meaning. In 
this lecture I will adopt the second point of view and try to explain what I am talking 
about. If you adopt the first point of view I leave you with the warning issued by Humpty- 
Dumpty 0. Words mean what I want them to mean, no more-no less. If they do a good 
job they come in on Saturday evening to collect their wages. 

It is only fair to close these general remarks with the candid acknowledgement that 
there is not universal agreement about what I will say. As an example it was stated in print 
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as recently as 1989 that the connection between (|1.1| ) and ( [Op was "not understood" ||. 
From the very title of these lectures it is clear that I do not share this view. 



3. Classical Statistical Mechanics 

We begin our considerations with classical statistical mechanics and in particular with 
a sketch of the theory of phase transitions and critical phenomena. These topics have been 
treated in various ways by other lecturers at this summer school so I will be concise. A 
primary intent is to use this procedure to make contact in sec. 5 with the quantum field 
theory of Q. 

The definition of the classical partition function Z and correlation functions in ( |1.3| ) is 
very general. However there are a variety of more special cases that are of great interest. In 
particular there is more than sufficient richness in the physics if we consider the following 
restrictions: 

1) The random variables <fi will be chosen to lie on the intersections of a lattice with 
L sites in each dimension (d = 1, 2, 3). 

2) The random variable (ft at each site r takes on a finite number of values N. 

3) The energy functional E({<p}) is translationally invariant with periodic boundary 
conditions. If pressed we will restrict the interaction between variables to be of finite range. 

Some examples of such energy functionals in two dimensions are: 

1) The nearest neighbor Ising model in a magnetic field 

L 

E = - ^ {E v (Tj ) kO'j+i,k + E h <jj^crj,k+i + Ha jjk } (3.1) 

where = ±1- 

2) The N state chiral Potts model 

L N-l 

E = - E {£>;,W+i,*r + E h n {a hk al k+1 Y + H n v% k } (3.2) 

j,k=l n=l 

where a^ k = 1. When iV = 2 we note that this reduces to (|3.1|) . 

The class of restrictions 1-3 is, of course, just one of many. The variables could lie 
on the bonds as well as the sites of the lattice. Continuous variables could be used and a 
familiar example of such a system is the n component classical Heisenberg magnet 

L 

E = - {vj,k ■ Vj + i t k + Vj >k ■ tfj,k+i} (3-3) 
j,k=l 
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where v 2 = 1. A continuum could be considered instead of a lattice and translational 
invariance can be broken. But for general discussion conditions 1-3 may be thought of for 
concreteness. 

We are interested in the thermodynamic limit when 

L — > oo, T fixed and positive. (3.4) 
In this limit we define the free energy per site in dimension d as 

/ = -kT lim T^lnZ (3.5) 

L— >oo L 

and we define the infinite volume correlations as 

lim < Oj > L (3.6) 



L- 



where < Oj >l is computed from (|1.3|) on the finite lattice. 

For a finite system the partition function Z of (|1.3| ) is an entire function of 1/T. 
However, once the infinite volume limit is taken the free energy (|3.5|) can have singularities. 
When the singularities occur for real positive T the system has phase transitions. 

These singularities in the thermodynamic limit are produced by the accumulation 
of zeroes of the finite size partition function as was first made clear in 1952 by Lee and 
Yang [|10| . In the general discussion it is necessary to distinguish between the cases where 



the zeroes pinch the real temperature axis at points and when entire segments of the axis 
are pinched. The pinching of an entire segment is known to happen in some cases where 
the interactions vary randomly from site to site and thus break translational invariance. 
We will consider some of these effects in later lectures. For the present, however, we will 
consider the case where the zeroes pinch the real temperature axis only at isolated points. 
This happens at what is called a critical point or a point of second order phase transition. 
If there are other external parameters in the problem, such as the magnetic field in the 
Ising model ( |3.1|) , the location of these points of phase transition can depend on these 
parameters. A plot of the location of the points of phase transition is called a phase 
diagram. One of the important tasks of statistical mechanics is the derivation of phase 
diagrams for realistic systems. 

Phase diagrams are quite dependent on the specific details of the system under con- 
sideration. Perhaps the most vivid illustration of this is the transition temperature of a 
superconductor. There is extreme interest in finding the material with the highest possible 
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transition temperature. While in principle this temperature is contained in the parti- 
tion function for the system under consideration, the practical truth is that theoretical 
statistical mechanics is far from being able to tell us how to make high temperature su- 
perconductors. 

However once the location of the points of phase transition are known there is a 
remarkable universality to the behavior of the system near the critical point. Thus although 
there are hundreds of superconductors there are only one or two types of superconductivity. 

Superconductivity is of course a quantum phenomenon and thus is not encompassed 
by the classical statistical mechanics we are considering here. But there are plenty of 
phase transitions in classical statistical mechanics as well. Familiar examples are the 
ferromagnetic transition to the magnetized state and the critical point in a liquid-gas phase 
diagram. Indeed, the discussion of []T(J on the close relation of these two phenomena is 
one of the foundations of the theory of phase transitions. Starting in 1944 with the exact 
work of Onsager on the Ising model JO]], extending through the 50's and into the early 60's 



with the development of various series expansions and culminating in the developments 
of the mid 60's with the invention of scaling theory, we have developed a very detailed 
picture of the physics of these classical critical phenomena. This is well expounded in the 
classic articles [0- [[nj. This picture is of fundamental importance for our discussion of 
the connection with quantum field theory. For concreteness we will develop this theory in 
the context of (|3.1|) and (|3.2[) . However the principles are generally valid. 

There are many phenomena which occur at an isolated critical temperature T c and the 
theory of critical phenomena relates them together. Some of the principal phenomena are 
a) the singularities in the free energy, b) the existence of spontaneous symmetry breaking, 
and c) the behavior of correlation functions at long distances. These are related through 
the construction of the scaling limit and scaling laws. We will discuss each of these points 
separately. 

3.1. Singularities in the Free Energy 

To discuss the singularities in the free energy it is convenient to define the specific 
heat c as 

d 2 f 

Then the simplest generic singularity the specific heat can have at a critical temperature 
T c is 

c~A a \T-T c \- a . (3.8) 
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The exponent a is referred to as a critical exponent. Somewhat more generally we recognize 
that the exponent a could be different if T — > T c from above or below. If the more general 
situation holds we let a be the exponent for T > T c and oi be the exponent for T < T c 
and write 

c~At(T-T c )- a (T^T+) and c ~ A~,(T C - T)~ a ' (T ^ T"). (3.9) 

3.2. Spontaneous Symmetry Breaking 

In the Ising model (371) the next simplest property to define after the specific heat 
at H = is the magnetization 

M(H) = lim < a , >l (3.10) 

L — >oo 

and for the iV state model similar order parameters are defined as 

M^(H) = lim < a™ > L . (3.11) 

L — >oo 

If H = the interaction ( |3.1[ ) is invariant under 

o~j,k —> ~°~j,k (3-12) 
and thus if M(H) is continuous at H = it follows that M(0) =0. If T > T c this is 
indeed the case. However, because we are considering the L — ► oo limit in the defini- 
tion ( |3.10|) or (|3.11| ), there is no reason that M(H) has to be continuous at H = and, 
indeed, for T < T c the continuity fails. Consequently for T < T c we define for ( |3.1| ) the 
spontaneous magnetization M_ (T) as 

M(T) = lim lim <a 00 >L ■ (3.13) 

Typically as T — > T c the spontaneous magnetization vanishes. Thus we define a second 
critical exponent (3 as 

M_(T) ~ A p {T c -TY as T^T~. (3.14) 

Another quantity related to the magnetization is the magnetic susceptibility \ 

dM(H) , 

X = dH \h=o- (3-15) 
This susceptibility also has singular behavior at T = T c and we parametrize this in terms 
of the exponents 7 and 7' as 

X~A+(T-T c )-t (T-T+) and X ^ A~,{T C - T)~^' (T —*T~). (3.16) 

Finally we consider the magnetization as a function of H at T = T c . This also has a 
singular behavior as H ~ which is parametrized as 

M(H, T = T c ) ~ A.sif 1 / 5 . (3.17) 

For the more general iV state model (|3.2|) the various exponents will depend on the 



index n of the order parameters (|3.11 
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3. 3. Correlations 

Not only do the bulk thermal properties of the system have singularities at T c but 
there are also corresponding phenomena in the correlation functions as well. Consider, for 



example, the two point correlation for (|3.1| ) in the infinite volume limit 

C(m,n) =< a ,oo-m,n > ■ (3.18) 

We write the coordinates as 

m = R sin 9, n = Rcos9. (3.19) 

When T < T c we find that as R — > oo the correlation approaches the limiting value of M_ 
exponentially as 

C(m, n) ~ M 2 _ (1 + °~^ T) e- R ^~ (e ' T) + . . .) (3.20) 
where £ _ (#,T) is called the correlation length. Similarly when T > T c 

C( m , n) ~ C+ ( g ' r ) e -fl/e + (g.T) + . . . . ( 3.2i) 

The correlation lengths depend on T and diverge as T — > T c , and thus we define exponents 
v and v' as 

+ T) ~ A U (9)(T — T c )~ v (T~^T+) and 

(3.22) 

C(e ) T)~A u ,(d)(T c -T)-' / (T^T+). 

The divergence of the correlation length as T — > T c is a signal that the forms ( |3.20| ) and 
( |3.21|) break down at T c and instead it is found that for T = T C the correlations decay as 
a power law 

C{m,n)~A^- (l^oc) (3.23) 
where d is the dimensionality of the system and r\ is called the anomalous dimension. 

3.4- Scaling Limit and Scaling Functions 

Of all the phenomena discussed above that happen at an isolated critical temperature 
T c the most important is the divergence of the correlation length £. The physical meaning of 
this divergence is that even though the physics is defined by a nearest neighbor interaction 
(on the scale that 1 is defined as the inter atomic spacing) the physical phenomena occur 
on a length scale £ which is much larger than this scale of definition. At the critical 
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temperature the physical scale is infinitely large compared to the scale of definition of 
the problem and thus it is most natural to renormalize our length scale from the atomic 
length to the observed physical length. Moreover on the atomic length scale the correlation 
lengths have an angular dependence on 9 which reflects the existence of the lattice on which 
the system was originally defined. However for many systems as T — > T c the angular 
dependence approaches that of an ellipse which can be made rotationally invariant by 
simply scaling the vertical and horizontal lengths with a different factor. Thus it is most 
desirable to define new renormalized lengths as 

x = n/f (0, T) and y = m/£(ic/2, T) (3.24) 

and take the limit 

T —> T c , n — > oo, m — > oo, with x and y fixed. (3.25) 

This particular limit is called mass renormalization and is part of what is called the scaling 
limit. 

If the limit (|3.25| ) were the only step involved the correlation ( |3.20| ) would vanish 
because the factor of the spontaneous magnetization vanishes and the factor C~(9,T) is 
found to go to a constant independent of T and 9. Consequently we also divide by M 2 
and define the renormalized Greens function as 

G(r)= lira MZ 2 C{m,n) (3.26) 

scaling 

where by lim sca i ing we mean ( |3.25| ) and r 2 = x 2 + y 2 . The process of dividing C(m, n) by 
M_ is called wave function renormalization. 

When T > T c there is a factor called M? similar to M_ which ( |3.21| ) may be divided 
by in order that the scaling limit exists. More generally we may consider correlations of 
n cr's. Thus we formulate the extension of the definition ( |3.26| ) to the general case of the 
scaling functions of the theory 

G(r u ...,r n )±= lim M± n < > (3.27) 

scaling 

where the subscript ± makes explicit that there are two different field theories to be 
constructed, one from above and one from below T c . 
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3.5. Scaling Laws 

Thus far the theory discussed may be considered as descriptive and all the exponents 
and functions introduced may be considered as independent subject only to the general 
requirements of thermodynamic stability. However if we make one additional assumption 
we find that this theory makes predictions about the relation of critical exponents to each 
other. These relations are known as scaling laws. 

The additional requirement is the assumption that there are no other length scales in 
the problem other that the atomic length scale of definition and the physical length scale 
of the correlation length and that these two scales join together smoothly. 

As a first example of the relations between exponents which this new assumption leads 
to we consider the relation between the scaling function ( |3.26| ) for T < T c and the lattice 
correlation ( |3.23| ) at T = T c . 

Considering the isotropic case for convenience and using the assumption of one length 
scale we extend ( |3. 19|) and ( |3.22|) as 



r = R{T c -TY /A v ,. (3.28) 
Then using again the assumption of one length scale we write ( |3.14| ) and ( |3.26| ) for T ~ T c 



as 

C(m, n) ~ A\{T C - T) 2l3 G(R(T c - T)"' /A v >). (3.29) 

The assumption of one length scale further says for large, but fixed R, that as T — > T c this 
form of G(m, n) must agree with the T = T c behavior ( |3.23| ). Thus by comparing these 
two expressions we find not only that it is required for the R dependence to match that 
for r ~ 



Go 

where Go is a constant but furthermore in order for the dependence on T c — T to match 
we need 

2(3 = v'(d-2 + r]). (3.31) 

This relation between critical exponents is called a scaling law. 

We may also use the one length scale assumption to find a relation with the exponent 
7' and an expression for A~, in terms of G(r). From ( |3.1(J| ) and ( |3.15| ) we may write 

*( T ) = £f E(< °o,ov m ,n > -Ml). (3.32) 
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In order to study the divergence in x as T — > T c it is sufficient to consider only the large 



R contributions to (|3.32|) . In this region the sum is approximated by an integral as 



/oo poo 
A vl {T c -TY v 'dx \ A V ,(T C -Ty v 'dy 
> . i -oo J — oo 

(3.33) 

/•2tt poo 

= A 2 u ,(T c -T)- 2v ' / dd rdr. 



o Jo 



Thus using ( |3.26[ ) we find that as T — > T c — 

p oo 

X (T) ~ (kT c )~ 1 A 2 ll (T c - Tf^-^Apn / r(G(r) - l)dr (3.34) 

Jo 

and hence comparison with ( |3.16| ) gives (written more generally in d dimensions) 

i = dv' - 2(3 (3.35) 

and (for d=2) 



poo 

AZ = (kT c )~ 1 A 2 jl A 2 3 2ir / r(G(r) - l)d 

Jo 



r. (3.36) 



This expression for 7' is another example of a scaling law. 

There are several other scaling law relations between exponents which are all obtained 
following the same line of reasoning using an assumption of one length scale. We refer the 



reader to |[L2|l - ||14jl for the details and here quote the final results valid for any dimension 

a' + 2(3 + 7' = 2, (3.37) 

dv' = 2- a (3.38) 

and 

The corresponding scaling laws for T > T c are obtained by replacing all primed exponents 
with unprimed exponents and replacing (3 with the appropriate exponent needed for the 
existence of the scaling limit ( |3.27| ) for T > T c . 
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3.6. Universality 

The final concept in the theory of critical phenomena that must be discussed is the 
concept of universality. In the examples of the Ising model and the chiral Potts model the 
defining formulae (|3.1|) and ( |3.2|) were written down with nearest neighbor interactions. 
But nothing in the subsequent discussion appears to depend on this assumption. Indeed it 
is thought that at the very least the interactions can be generalized to an arbitrary number 
of finite range (ferromagnetic) interactions and the only change in the above discussion will 
be that the amplitude factors and T c become functions of all of the interaction constants. 
However the exponents a, /3, 7, 5, 77, v and the corresponding primed exponents are inde- 
pendent of all of these interaction constants. More generally all of the scaling functions are 
independent of the interaction constants. This independence is referred to as universality. 

3.7. Discussion 

This general picture of the scaling theory of the critical point was developed over 
25 years ago. At the time it was originally put forward the proponents were exceedingly 
cautious. The nature of the one length scale assumption was clearly investigated and the 
consequences of singularities more complicated than the pure power laws were examined. 
But a quarter century is a long time and it is appropriate to close this section by surveying 
the subsequent progress made in verifying this theory. 

The first very significant piece of information is that in the last 25 years it has been 
verified that the two dimensional Ising model at H = totally confirms all aspects of the 
theory. Thus the class of scaling theories is definitely not empty. The full explanation of 
these results will be given in lecture two. 

There is no other theory as well studied as the two dimensional Ising model and 
in three dimensions there are only numerical confirmations of the scaling laws available. 
However, the current status here is that after 25 years of most serious endeavor no known 
counter example has been found. 

In two dimensions there exist a plethora of models for which partial results on the 
scaling behavior are known. In all known cases the simple power law forms hold (unless 
the exponent is an integer in which case a logarithm may appear). Moreover whenever 
sufficient exponents have been computed to test a scaling law the prediction has always 
been verified. However for no model other than the Ising model have the scaling functions 
been fully characterized. 
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The verification of the concept of universality is less studied. Even for the Ising model 
exact computations have never been done if the interaction bonds do not form a planar 
graph. However there is absolutely no evidence that universality fails. This is one of those 
cases in physics where it would be most desirable to have a precise mathematical theorem. 



4. Quantum Statistical Mechanics 

We next turn to quantum statistical mechanics as defined by (|1.6|) . As with classical 
statistical mechanics it helps to consider specific systems in order to illustrate the problems 
involved. An example that will possess more that sufficient richness is the XYZ chain in 
an external field 

h = - ^(./'.s;-;.s;; . , + .r<s ! ;s ! <. , + J z s*s* +1 + hs z 3 ) (4.1) 

where S l can be either the three Pauli spin matrices a 1 or the spin S representation of 
the rotation group. When J x = J y = J z the resulting chain is called the spin S isotropic 
Heisenberg magnet. 

As in the case of classical statistical mechanics, phase transitions exist in quantum 
systems. Unfortunately our knowledge of the details of these transitions is much less 
than that of their classical counterparts. Thus although the theory of critical phenomena 
developed in the previous section is just as applicable to the quantum as to the classical case 
and even though there is series expansion data for critical exponents of several quantum 
magnets in three dimensions we will not pursue the topic further in this lecture. 

However there is one new feature in quantum statistical mechanics which is very im- 
portant to consider. Namely the quantum dynamics which is obtained from time dependent 
correlation functions such as 

C J (n,t;T) = Tr(e- ltH S J Q e ltH Sle- H / kT )/Tre- H / kT (4.2) 

Just as the theory of equilibrium critical phenomena starts from the behavior of the 
correlations for large spatial separations, so the theory of transport properties starts from 
the behavior of these time dependent correlations for large times. This theory originates 
in the work of Bloembergen [15|] and de Gennes [16|. We here follow the treatment of 



TTJ and . 
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It will make a great deal of difference in the asymptotic behavior whether or not 
there is a conservation law. We will concentrate here on the case where there is such a 
conservation law of the form 

oo 

C*(n,t) = l. (4.3) 



n= — oo 



Such a conservation law holds for the XYZ model ( J4.1[ ) if J x = J y . For this case we define 
the spatial variance 

oo 

a 2 (t) = £ C z (n,t)n 2 (4.4) 

n= — oo 

and the spatial Fourier transform 

oo 

I(k,t)= C z (n,t)cos(kn). (4.5) 

n= — oo 

The process of spin diffusion is said to occur if for small k and large t we have 

J(M) ~ Ae- DkH (4.6) 

where < D < oo. This asymptotic behavior implies the following long time behaviors 
for the autocorrelation function 

/*" dk 
—I(k,t) ~ (AnDt)- 1 / 2 (4.7) 

and the spatial variance 



° 2 (t) = -^f l U=o~2£*. (U 



d 2 I(k, t) 

These long time behaviors play the same role as does the exponential decay of the 
spatial correlations discussed in the previous section and as before we need to know the 



validity of (|4.6|). In 1976 it was stated [17 that "this results rather naturally from the 



phenomenological theories of irreversible processes, and can be more or less proved by 



modern statistical mechanics, for instance by Mori's methods ||19|1 ". However, despite the 
optimism of this claim the fact is that this spin diffusion form has never been demonstrated 
for any system and the only exact computations which do exist have a different asymptotic 
behavior altogether. This is not to say that the spin diffusion form does not occur, but 
it does mean that our knowledge of the subject is not nearly refined as our knowledge 
of critical phenomena in classical systems. Consequently, even though these quantum 
time dependent correlation functions are naturally connected with finite temperature field 
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theories, we will not pursue the connection further in this lecture. We will consider these 
phenomena at length in lecture four. 



5. Quantum Field Theory 



The formula for thermal averages (|1.3| ), which was the basis for sec. 3, was found 
by Gibbs and Boltzmann in the 19th century. It forms the first chapter of any course 
on statistical mechanics and all our treatment of thermal phenomena follows from ( |1.3| ). 
There is never any argument about the starting point of equilibrium statistical mechanics. 

On the other hand the formula ( |1 . 1| ) has no such universal status in quantum field 
theory. Consequently it may be useful to discuss where ( |1 . 1|) comes from. 

It is here, surely, that the warning of Humpty-Dumpty must be kept in mind because 
the words quantum field theory mean many things to many people. Some of the different 
usages are exemplified by the following list: 

1) Canonical QFT 2) Lagrangian QFT 

3) Relativistic QFT 4) Perturbative QFT 

5) Renormalized QFT 6) Finite Temperature QFT 

8) Topological QFT 
10) Gauge QFT 
12) S Matrix Bootstrap QFT 
13) Path Integral QFT 
The formula (1.1) may be taken as the definition of the last item on the list, but there 
are literally tens of thousands of papers written about the other twelve items and this list 
is far from complete. In fact, there are so many different usages of the words "quantum 
field theory" that it often useful to regard the words QFT by themselves as not specifying 
anything at all until the qualifying adjective is added. This warning saves untold hours of 
discussion which communicates nothing because the participants are all talking at cross 
purposes. It is certainly not the case that all books on quantum field theory start with 



7) Conformal QFT 
9) Lattice QFT 
11) Axiomatic QFT 



[1.1) as the definition. Some books |2(J [^TJ put (1.1) in the very title, some make use of 



it in the text in places ranging from the very beginning to near the end H22| - |2ql and some 
do not use it at all |26] p7 [. This lack of a common starting point marks a fundamental 
distinction between statistical mechanics and quantum field theory. 

With this being said, however, it is certainly true that the first twelve items are not 
totally unrelated to the thirteenth. Therefore I will begin this section by sketching the 
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history of the development of formula (1.1). As in the previous two sections it is useful to 
have several examples in mind to make the considerations concrete. Three such examples 
of Lagrangians are |23[ 



1) The n vector model 



d n ? 9 



k=i k 

This model can be viewed as a generalization of the n component Heisenberg model 
(which in quantum field theory is often called the non-linear sigma model). 
2) Scalar Electrodynamics 



4 

j,k=l 



d 

+ 5}(0J + ieAjW^idj + ieAj)<j>] + mV</> + f 

3=1 



(5.2) 



3) Quantum Electrodynamics 

L = \ E ^ Ak ~ dkA rf + ^{d-Af + ^{d- 1 -eiA- 1 - m)V (5.3) 

j,k=l 

where 7j are the Dirac gamma matrices. 

Quantum field theory did not originate with the discovery of (1.1). Instead quantum 
field theory developed from classical field theory by taking the equations of motion of a 
classical field (of which the Maxwell's equations coupled to electrons was by far and away 
the most important application) and replacing classical fields by quantum operators which 
obey commutation or anti-commutation relations. From these field equations Green's 
functions were computed and from this such triumphs as the anomalous magnetic moment 
of the electron came. 

Classical field theory can be obtained from a classical action and as early as the work 
of Dirac in 1933 [28] and the work of Feynman [29] in 1948 it was understood how to 
extend this to quantum mechanics using a path integral of the Minkowski form ( |1.4| ) . The 
action principle was extended from quantum mechanics to quantum field theory shortly 
thereafter by Schwinger |30|| . At almost the same time it was realized by Dyson |[31|| , and 
later by Wick |3!| , that great advantages could be obtained if one considered the analytic 
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continuation to Euclidean space and the Euclidean action approach of (1.1) was considered 



by Schwinger |33[ in the late 50's. 

It may now be asked what is the origin of the action principle. This can be answered 
by first forming the partition function Z of (1.1). If h is considered to go to zero and if we 
choose to evaluate the functional integral for Z by the method of steepest descents where 
the action is expanded about classical configurations which satisfy 

SL 



8(ft(x) 



(5.4) 



then the results of Schwinger's action principle and of canonical quantization are the same. 
If in addition external sources are considered then the partition function can be regarded 
as the generating function for the Greens functions of the quantum theory. In this con- 
struction h sets the scale of the commutation relations and this scale is identified with the 
experimentally determined Planck's constant. 

Thus far the path integral (1.1) has been used as a device for doing computations 
which could and indeed were first performed by other methods. 

However in this procedure we have been rather sloppy because questions of existence 
and convergence have been ignored. In order to make direct sense of (1.1) some sort of 
additional definitions need to be made. 

For the scalar field theory 1) it is sufficient to replace the fields (ft(x) defined on 
continuum space by fields defined only on the lattice sites and to replace derivatives by 
differences which amounts to the replacement 

(d k (ft) 2 - -0(i)0(xO (5.5) 

where x and x' are nearest neighbors. 

For the gauge field theories 2) and 3) the lattice definition was first given in |34|| . The 
vector potential %A is replaced by a compact group element such as e which is defined 
on the links of the lattice and the matter fields are (ft or ift are defined on the sites of the 
lattice. Then, calling the gauge group element on the link between sites j and k Uj ik the 
Lagrangian 2) is replaced by 

^-TrU jik U k)l U lim U mjj + e^(x j )U jjk (ft(x k ) + j(<f>Uf (5-6) 

Zg Q Z 4 

where in the first term the sites j, k, /, m lie on the vertices of a square (plaquette). 
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But now that we have reproduced the statistical mechanical starting point we are 
forced, if we want to regain a continuous result, to follow the procedure that leads to the 
scaling limit in statistical mechanics and find those points in the phase diagram, including 
the variable h, where the system undergoes a phase transition and at those points to 
construct a scaling limit. If the point h = is a point of phase transition then we will regain 
the situation we started from with canonical quantization. But in statistical mechanics 
there are a host of physically interesting systems where the critical temperature is not at 
zero. The corresponding field theories will be outside the realm of those considered by 
canonical quantization. 

We may now examine the concept of perturbative QFT as derived from (1.1). In 
principle one obtains a field theory from (1.1) by finding the points of phase transition and 
constructing a scaling theory in the vicinity. In perturbative QFT the action is written 
as Sq + XSi where Sq has a phase transition at h = 0. An expansion in A is then made 
assuming that the full theory still has a phase transition at h = 0. Actions for which 
this procedure works are called renormalizible and the resulting expansion is called a 
renormalized perturbation expansion. These renormalized perturbation expansions form 
a large part of the literature on quantum field theory. 

In these renormalized perturbation expansions the physical content of particles of the 
full theory is the same as the content of the "free" theory So (up to possible questions 
of bound states). The particle spectrum of Sq is called bare as opposed to the genuine 
physical particle content of the full S which are called dressed or renormalized in particle 
physics. 



But starting in the mid 60's, certainly with the seminal paper of Gell-Mann |p5[ , it 
has been realized that for the gauge theory of strong interactions (QCD) that neither the 
gluons (gauge bosons) nor the quarks (bare fermions) which are the fields in S actually 
show up as physical particles of the full theory. The word for this is confinement and 
the lattice gauge theory of |34j] is one way of addressing the problem. (Any reader who 
has not in fact read |?5[ is strongly urged to read the last sentence and to marvel at the 
subtlety of the English language. Even Humpty-Dumpty must stand in awe of an author 
who can introduce into a physics paper a word whose origin is explained by a footnote to 
Finnegan's Wake and then suggest an experiment to determine if these objects are real 
and if they exist.) 

The need to consider field theories with confinement is a most compelling reason to 
extend the notion of QFT beyond the perturbation theory fixed point at h = 0. 
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Consequently there has been a movement in the past 30 years to change our point of 
view and to take (1.1) as the starting point for QFT and to allow for all possible points of 
phase transition, not just h = 0. 

This change of starting point for QFT now makes QFT and SM into identical subjects. 
This is natural for a statistical mechanician but from the historical perspective of QFT we 
must accept the fact that if the point of phase transition is not h = 0, there will no longer 
be a classical limit of the QFT. Such a field theory may rightly be called fully quantum 
mechanical. 

In more old fashioned terms we are accepting the fact that for field theories with 
points of phase transition not at h = there is no correspondence principle. 

In the past 20 years we have had great success in two dimensions in finding, construct- 
ing and in many cases solving these quantum field theories with no classical limit. Indeed, 
most of the conformal and topological field theories discussed by other lecturers at this 
summer school have this property. 

Moreover once it is understood that fully quantum mechanical quantum field theories 
exist and that at least in two dimensions some of them can be explicitly seen to have 
confined degrees of freedom, it is not possible to accept as dogma the predictions about 
high energy scattering made solely on the basis of renormalized perturbation theory. The 
development of QFT from quark confinement to string theory is the attempt to understand 
physics beyond the perturbation theory of the 60's. 

At this point it is most helpful if a very candid admission is made. If one reads the 
papers on the path integral foundation of quantum field theory just discussed above |28 - 



||33|| it is very clear that this procedure of looking for points of phase transition in the 
Euclidean path integral is NOT what these authors had in mind. We have quite deliber- 
ately extended the meaning of the words quantum field theory beyond what the founders 
of quantum electrodynamics had in mind by equating quantum field theory with phase 
transitions in statistical mechanics. The reader has every reason to ask who it was that 
first made this new extension. 

It is one of the curious mysteries of physics that there does not seem to be any one 
person or paper to whom this idea can be unambiguously attributed. In his 1974 paper 



Wilson []33[] writes as if it were well known and refers the reader to the paper of Wilson 
and Kogut |3(J. But in |3(J the authors also do not claim that they are inventing a new 
principle. The statement of many of the ideas relating field theory and critical phenomena 
in a perturbative setting is given several years earlier [3~71]- HU where the focus seems not 
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to be on the invention of a new principle but on rather heavy diagrammatic computations 
and the relation with field theory is treated as rather obvious. Indeed in [^(J the Ising 
model is discussed and the Ising model had been discussed in terms of a field theory of 



fermions as early as 1949 [g4[] . 

My conclusion from this is that the emergence of the identity of QFT as defined by 
(1.1) and statistical mechanics as defined by ( |1.3| ) emerged from the community as a 
whole as something that once you saw what was going on you truly decided that it was 
so obvious and basic that it was implied in all that had gone before. It is perhaps fitting 
that a thermodynamic principle has a collective origin. 

However, the lack of a definitive statement of origin leads to confusion even to the 
present day. If you ask many physicists if there can be a quantum field theory without 
a corresponding classical limit you will often get a blank stare of incomprehension. And 
yet it tortures the brain to force the Ising model ( |3.1[ ) into the language of classical field 
theory and a stationary action principle. In particular, as explained in sec. 3 there are two 
field theories that can be constructed near the T c of the Ising model: one from above T c 
and one from below T c . Even worse is the fact that there are actually an infinite number 
of field theories that can be constructed as H — > and T — > T c depending on the value of 
h = H/\T -T c \ 15 / 8 . One way in which it has been attempted to preserve the notion of a 
classical limit is to consider the Ising model as an n = 1 non linear sigma model (|5.1| ) where 
the interaction term is A(</> 2 — l) 2 in the limit where A — > oo and to develop a perturbation 
scheme in which e = 4 — d is a small parameter. This was originally suggested in [^5j and 
is reviewed in detail in [l|]. However the e expansion does not converge and thus far 
it has not been possible to use it to regain the exact results known in two dimensions. 
Consequently I will here adopt the view that it is far better to abandon the notion of a 
correspondence principle rather than to turn an exactly solvable model into an insolvable 
mess merely to preserve a fiction that indeed may not represent all the physics. 

It finally remains to consider how to obtain results in the physical Minkowski metric 
from the Euclidean path integral (1.1). Formally this is done by analytically continuing 
one of the Euclidean variables, say xq, to it. The mathematical questions that need to be 
addressed are 1) the circumstances under which this continuation is possible and 2) what 
the relation is of the resulting continuation with correlations computed directly from the 
Minkowski path integral ( |1.4j ). Fortunately these questions have been extensively investi- 
gated in axiomatic quantum field theory and are well discussed in the literature [pTf , [[I6| . 
Consequently we can be brief. 
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For a Euclidean QFT to be acceptable it is not enough that the limits exist but 
that the limiting correlations satisfy several properties called the axioms of Osterwalder 



and Schrader [46|. In non technical language these axioms are: EO) Temperedness; El) 
Rotational invariance; E2) Reflection (Osterwalder-Schrader) positivity; E3) Symmetry; 
E4) Cluster property; 

The virtue of the OS axioms is that they hold for the Euclidean QFT if and only if 
the following Wightman axioms ||47|| - ||49|1 hold for the Minkowski correlation functions: 
MO) Temperedness; Ml) Lorentz Invariance; M2) Positivity; M3) Local commutativity; 
M4) Cluster property; M5) The existence of a lower bound on the energy spectrum. 

For further discussion of these axioms (especially for the somewhat technical expla- 



nation of temperedness) see ref [f46| and ||21|1 . 

The verification that the correlations of a Euclidean QFT satisfy the OS axioms can be 
taken as guaranteeing that the analytic continuation from Euclidean to Minkowski space is 
possible. Thus for our purposes we will consider that whenever they are both defined that 
(1.1) and ( |1.4| ) are equivalent. It may, however, be that there are cases where ( |1.4| ) is 
not necessarily well defined. In that case we will use (1.1) and analytic continuation as the 
definition of what is meant by a Minkowski QFT. 



6. The Need to Exist 

One of the enduring contributions made by Kant in The Critique of Pure Reason is the 
demonstration that the mere making of a definition does not by itself ensure the existence 
of the thing defined. Kant raises this most important distinction in his discussion of the 
ontological proof given by Descartes in the Discourse on Method for the existence of God. 
However the truth of this distinction is universal and since Kant's time the existence proof 
has become a very important part of mathematics. 

Consequently the definitions of path integral quantum field theory given above will 
only be useful if we can demonstrate that the category of such field theories is not empty. 

In statistical mechanics a great deal of attention has been given to the determination 
of restrictions which must be placed on the classical energy functional E((p) or the quantum 
Hamiltonian H in order for the thermodynamic limit of the correlation functions ( |1.3| ) or 
( |1.6| ) to exist [5D]-[5S]. To prove existence it is necessary that the interactions are not 



so attractive at short distances that the system collapses and that at long distances the 
interactions are not so repulsive that the system explodes. Somewhat more precisely the 
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requirement of no short distance collapse will be ensured for a quantum system if the 
ground state energy has the lower bound of 



E GS > -AN (6.1) 

where iV is the number of particles in the system and A is a positive constant. 

It is one of the triumphs of mathematical physics that the lower bound ( |6.1| ) has been 
proven for the system of electrons and protons interacting via the non-relativistic Coulomb 



interaction ||52|| . ||54]| and ||55|| . For the existence of the bound it is essential for the electrons 
and protons to obey Fermi statistics. A system consisting of only charged bosons is known 
to have an upper bound on the ground state energy of |HJ [^2| 



E GS < ~A X N 7 ^ (6.2) 

which manifestly violates (|6.1| ). Such a system is unstable against collapse and the laws of 
statistical mechanics (|1.3|) are not directly useful in its description. 

There is an identical need for the demonstration of existence for quantum field theory 
as defined by (1.1). However similar detailed existence proofs have not been carried out 
for all actions commonly believed to be interesting. But just as in the case of statistical 
mechanics the need for existence places restrictions on the actions which can be considered 
and in particular a lower bound on the action similar to (|6.1| ) is required. 

If it were the case that all the actions which are considered interesting in quantum 
field theory had such lower bounds then it would not be particularly appropriate to dwell 
on this point. However, this is the not case. The most famous case in point is the attempt 
to quantize the Euclidean version of the Einstein action of general relativity 

Se = -ltd j ' Rgl,2dlx < 6 - 3 > 

where R is the curvature and g is the determinant of the metric tensor g^^. This classical 
action is indeed not bounded below and thus the definition of QFT given above based on 



(1.1) will not apply. There have been attempts [56]— [58] to introduce extra restrictions 
into the Euclidean path integral to exclude the regions of path space which make (1.1) 
diverge. There have also been proposals |59|] for different Euclidean actions obtained from 



a more sophisticated analytic continuation procedure than used in non gravitational field 



theories. Some of these quantization procedures seem to lead [^0[ to further difficulties 



Moreover, even if such a prescription could be found there is a further problem of how 
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to pass from the Euclidean to the Minkowski metric. In a general relativistic theory the 
correlations are certainly not expected be to be rotationally or Lorentz invariant and thus 
there is no known equivalent of the Osterwalder-Schrader axioms [ |6| . Consequently there 
is no reason to believe that there is enough analyticity to continue from a Euclidean to a 
Minkowski region. Indeed the global question of what is meant by Minkowski is obscure 
for quantum gravity [pD|. 

Of course it is formally possible to construct a perturbation expansion about the point 
h = if (1.1) is used in an uncritical fashion and the classical solution is asymptotically 
flat |K| |nj , but the resulting perturbation theory is readily shown to be nonrenormalizible 
and is thus not particularly suitable to be taken as a definition. 

One physical effect that lies at the bottom of the difficulty of using (1.1) to provide 
a theory of pure quantum gravity is that gravity is a purely attractive bosonic theory and 
this attraction at short distances leads to precisely the sort of collapse that occurred in 
the system of nonrelativistic charged bosons. 

Because gravity is always with us and has been studied for literally 3000 years it is no 
small criticism that this interaction is not well incorporated into the definition (1.1). I am 
also well aware that there is a school of thought that argues that the only way to obtain a 
consistent theory of elementary particles is to find a theory which does incorporate gravity. 
Many attempts have made in addition to []56| - f)9| to either 1) modify the classical action 
to make the path integral exist, 2) add sufficient fermionic fields (usually of spin 3/2) to 
provide enough repulsion to produce existence of the path integral or 3) abandon (1.1) 
altogether and to develop new principles. Some of these attempts are summarized in 
||60|| but it is fair to say that none has been successful on all counts. 

It is thus an open problem to find a definition of QFT that will incorporate gravity. 
But following the philosophy of Humpty-Dumpty I will not let this stop me from defining 
what / mean by quantum field theory. Consequently in these lectures the words quantum 
field theory will always mean the limiting theory constructed from a lattice regularized 
path integral (1.1). No more-No less. 



7. Particles 

I have now given some sort of meaning to the formal connection between quantum 
field theory and statistical mechanics based on ( |1.1| ) and ( |1.3| ) and I have sketched many 
of the physical phenomena measured in the statistical mechanical part of the connection. 
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But in terms of quantum field theory I have said nothing at all about the experiments 
that quantum field theory was invented to study in the first place; namely the study of 
the properties of particles. 

The word particle is one that has many connotations. Particles are conceived of as 
having properties such as momentum, charge, mass, spin and a variety of other quantum 
numbers given in the particle data tables. There is the distinction between stable and 
unstable particles and there is the question of how unstable a particle can be and still be 
dignified by that name. Particle may be "elementary" such as neutrinos and electrons or 
"composite" such as nuclei and we may argue over which of the two categories neutrons 
and protons are supposed to be in. Particles such as neutrons and protons have some 
sort of intrinsic size or extension. To complete the relation of quantum field theory and 
statistical mechanics we must say something about the statistical mechanical counterparts 
of these objects. 

In this lecture I will not address all of these different properties connected with the 
word particle. In particular I will not consider what should be meant either by the word 
unstable or the word extended. Thus the discussion of this section will only begin to study 
this most fundamental of concepts. 

I will begin with the concept of a particle like energy spectrum of a Hamiltonian. Such 
a spectrum has the property that in infinite space all the excited state energies are given 
by 

S m a 

E ex -E GS = J2 E e ^ t 7 - 1 ) 

a=l j a = l, rules 

and all the momenta of the states are given by 

S m a 

= £ E ( 7 - 2 ) 

a = l j a =l, rules 

where S is the number of different species of particles, m a is the number of particles of 
type a in the state, e a (P) is the single particle energy of the particle of type a, and the 
subscript rules indicates the rules with which the excitations can be combined. One such 
example of the combination rules is 

Pj a + Pl a if 3* ^ la- (7.3) 
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This is a called a Pauli exclusion rule and particles which satisfy it are called fermions and 
said to satisfy Fermi statistics. Another example of such combination a rule is 



m 1 +m 2 = (mod 2) (7.4) 

which expresses the fact that particles of type 1 and 2 must be created in pairs. When this 
type of rule exists there is no one particle state for either type of particle and the single 
particle energy must be inferred out of the observed multiparticle energy levels. 

In a relativistic quantum field theory the single particle energies are of the form 

e a (p) = (M*+f) 1 / 2 (7.5) 

and the statement that all the states of the system satisfy ( |7.1| ) is referred to as asymptotic 
completeness |49] . It is an open question as to whether the Wightman axioms are sufficient 



to guarantee this form. 

In statistical mechanics this type of spectrum is referred to as a quasi-particle spectrum 
The prefix "quasi" being attached to indicate the feeling that somehow magnons and 
phonons are different from protons and electrons and that single particle energy levels 
other than (|7.5| ) are allowed. 

This particle spectrum has a major impact on the Greens functions. This is best 
described by considering the Fourier transform of the correlation functions (|1 . 1[) or ( |1.4|) . 
For concreteness consider the Fourier transform of a Euclidean two point function 

G(k) = J dx C{xy % ' 3 (7.6) 

where C{x) is obtained from the Minkowski correlation function by Wick rotation. Then 
from (|7.1| ) and ( ff75|) it follows that G{k) depends only on k 2 and has poles at 

P = -Ml (7.7) 

and has multi-particle threshold singularities at all those combinations 

^ = "($>aM a ) 2 (7.8) 

a 

where n a takes on all integer values consistent with the rules of combination. 
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We will take the spectral condition and the corresponding statement on the lo- 

cation of singularities in the Greens functions (|7.7|) - (|7.8|) as the definition of the existence 
of particles in a relativistic quantum field theory. 

To explain further the relation of particles to Greens functions in statistical mechanics 
we introduce the concept of the transfer matrix T (considered in two dimensions for 
concreteness) . Consider a system with nearest neighbor interactions on a square lattice 
where the energy is of the form 

1=1 3=1 

where E v ((frjj, <j>j+ij) (E {<j>j,i,<fij,i+i)) are the vertical (horizontal) interactions. Then 
letting {</>} stand for the collection of variables 0^ in row j and defining 

rph r - V ih E h ((j}i,(j)i + 1 )/kT f>-,-\c\\ 

and 

a row to row transfer matrix may be defined either as 

t r _ T v T h {7.12) 

or more symmetrically as 

rpR 1 /rph^l/2rpv ^g-j 

The partition function (|1.3|) with periodic boundary conditions in both the vertical and 
horizontal directions is then written in terms of either of these transfer matrices T as 

Z = TrT L v. (7.14) 

It is also useful to consider a diagonal to diagonal transfer matrix 

T W,{4>'} ^e^^e^W/^ (7.15) 
l=l 

Then if the boundary conditions on the lattice are suitably shifted the partition function 
is still given by ( |7.14j ). 
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In general transfer matrices are diagonalizable. Thus calling A m the eigenvalues of T 
we may write ( f7.14j) as 

z = J2 x ™ ( 7 - 16 ) 

m 

and in the thermodynamic limit the free energy ( |3.5| ) is written in terms of the maximum 
eigenvalue Ao of T as 

f = -kT lim — lnA . (7.17) 

The correlation functions may also be expressed in terms of the transfer matrix. As 
an example consider the two point correlation of the operator aji of the chiral Potts model 
By definition 

states 

which is written in terms of the transfer matrix as 

< ctVIm' >= ^TrafTJ'-iaPT^-^-^ (7.19) 

(where for definiteness we have chosen j' > j). In the basis where the spins take on 
the values e 27Tir / N with r = 0, 1, . . . , N — 1 the operator o\ is 

o- l = J®...® v 5^®...® (7.20) 

site I 

where S z is the N x N matrix with elements 5 r ye 2 ' K ' ir t N Thus, writing T in terms of its 
eigenvectors \m > and eigenvalues A m we find in the thermodynamic limit 

< vlivfv >= < °\ a ?\ m > (^m/Ao)°"" j) < m|a;, n |0 > . (7.21) 



Corresponding to the quasi-particle property (|7.1| ) of the energies there is a corre- 
sponding quasi-particle form for the eigenvalues of the transfer matrix 

S m a 

Aex/A =n n a «^)- ( ? - 22 ) 

a=l j a =l, rules 

Moreover, for a translationally invariant system the operator is obtained from the 
operator ao by a lattice translation operator whose eigenvalues are e lkPe:c where P ex is 
given by ( |7.2|) . Thus the correlation ( |7.21| ) is expressed as 

5 oo 

< «>?r >=EE E < °l <r o l{Pi«} >< feJkolO > 

a=lm Q =0{p jQ } (7.23) 

eELiE^Li^^'-Opja+U'-iJWAaCp^))] 
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This representation of the correlation function is the statistical mechanical analogue of the 
Kallen-Lehmann representation of quantum field theory ]23|] where lnX a (p) plays the role 
of the single particle energy levels. If the scaling limit exists and is relativistically invariant 
with S particles of masses M a we obtain the two point Greens function as 

G n (r) = lim M~ 2 < a>*,> > 

scaling 
S oo 

= E E < °kol{Pi.} >< iPjJWo n \0 > (7-24) 

a=l m a =o 

e Ef =1 =1 [ir * Pj <* +rv (M " +p l<* )1/2] 

Here the matrix elements < 0|o"q |{p Ja ,} > are referred to as form factors. In particular for 
the correlations to be rotationally invariant the single particle form factors < 0|(7q \pj a > 
must be constants. Then the Fourier transform ( |7.6| ) will have poles at ( [7.7| ) and we 
see that the two notions we have introduced for the masses M a of the particles coincide. 

But in applications of field theory to particle physics the notion of particle spectrum 
and masses of particles are not sufficient to describe the experiments. We must further 
be able to describe the scattering of particles. Particle scattering has been at the core of 
elementary particle physics experiments for most of this century and I will presume the 
reader has some acquaintance with the subject. For concreteness I will follow the treatment 



of IE3f 



The fundamental object in scattering theory |23| is the S matrix which represents 
the transition amplitude from an "in" state with n- m particles of momenta {p a } to an 
"out" state with n ou t particles of momenta {p' a }. One of the crucial results is that the S 
matrix elements, which are supposed to be what are measured in laboratory experiments, 
are expressed as the residues of Greens functions when the external momenta satisfy the 
mass shell condition ( |7.7| ). This result is the LSZ reduction theorem ||23|| . 



However, when the attempt is made to relate the LSZ reduction theorem to the 
discussion of Greens functions given above there are two problems which cannot be avoided. 

The first problem is that there may be many poles in a two point function and thus 
there does not need to be a one to one relation between particles and fields. This is an old 



observation [5^] and for such systems the notions of "creation operators" and "free fields" 
which are used in the derivation of the LSZ theorem need a more careful treatment. To 
make a "particle" to be observable we would like, in some theoretical sense at least, to be 
able to find an operator that produces only that single particle. The definition of mass 
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and spectrum made above do not provide an answer to this question. This phenomenon of 
many poles in a two point function occurs in the field theories constructed from the Ising 
model in the limit H — > 0, T — > T c with H/\T — T c | 15 / 8 held fixed and will be discussed in 
detail in the third lecture. 

The second problem is the question of which Greens functions are to be used in the 
LSZ reduction formula. In the iV state chiral Potts model explicitly discussed above there 
were N — l different two point functions which were considered, all of which are constructed 
from the same energy spectrum. The same is true for the four and higher point functions. 
It would thus appear that I can construct several different S matrices for the same particles 
by use of the LSZ reduction. This would seem not to be compatible with the notion of a 
particle as a thing in itself and consequently it can be argued that the consistency of the 
particle interpretation will force an identity between these several different constructions. 
However I know of no proof of this for any model and, in fact, do not even know where 
the point is discussed. 

If in fact there are different S matrices obtained from different Greens functions we 
would be forced to somewhat modify our interpretations. However it does seem to reflect 
reality to adopt the point of view that in any experiment we will only say that a particle is 
observed if it is detected in some way. The detection process must refer to the interaction 
with some classical external apparatus. The word "observable" is used indiscriminately in 
speaking of the operators Oj used in the correlation functions, but in actual practice we 
can only observe those operators for which apparatus can be manufactured that will couple 
to it. It would seem that at least in terms of Greens functions that it is an act of some 
philosophic subtlety to speak of a particle in the abstract without considering a concrete 
measurement and it does not seem excluded that different measuring devices could in fact 
detect different objects. 

Particles are not things of direct sensory perception but are rather philosophic ab- 
stractions. Their existence is the existence of the abstract and no matter how appealing 
it is at times to think of them as small billiard balls such a concept is not adequate. We 
will return to this in subsequent lectures. 



8. Statistical Mechanics and Quantum Field Theory Compared 

I have now provided one answer to the question of the connection of statistical mechan- 



ics and quantum field theory via the path integral construction (1.1) and (1.3). However 
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to make completely clear the relation between the two subjects I will close this lecture by 
pointing out that the scaling limit has a very different status in statistical mechanics than 
it does in quantum field theory. 

In the theory of magnets and critical fluids the notion of the atomic length scale is 
very real and physical. Moreover the critical temperature T c is directly observable. There 
are many measurements of specific heats, order parameters, and magnetic susceptibilities 
for arbitrary temperatures away from T c . These non critical (and nonuniversal) properties 
of real statistical mechanical systems are very familiar. What is difficult to measure in a 
magnet or a fluid are the properties in the critical region. Because critical theory is by def- 
inition a limiting theory the question of how close to T c you have to go before the formulae 
of the theory may be applied is of very great practical concern. The scaling functions are 
almost never measured in detail because of the difficulties in precisely controlling T — T c . 
For an in depth discussion of these problems the reader is referred to |)T . 



In quantum field theory, on the other hand the situation is exactly reversed. Here the 
scaling functions are the measured scattering amplitudes whereas the order parameters and 
the correlation lengths are the renormalization constants which can never be measured. 
Indeed because the bare theory and the physical theory are cut off by an infinite renormal- 
ization of scale it is not totally clear how unique the concept of the bare fields are. Even if 
in some sense it is unique it is most doubtful that it can ever be reconstructed from a finite 
number of measurements. We will return to this question of a practical reconstruction 
theory in the next lecture. 

The spirit of the relation of the scaling to the bare theory in both SM and QFT is 
captured by the old children's rhyme: 

I met a man upon the stairs. 

The little man who was not there. 

He wasn't there again today. 

My God, I wish he'd go away. 
The difference between the subjects is that in statistical mechanics the critical theory 
is the little man whereas in quantum field theory he is the bare theory. 
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Part II. Ising Field Theory for H = 



Abstract 



We discuss the theory of the two dimensional Ising model at zero magnetic field in 
both the homogeneous and inhomogeneous cases. For the homogeneous case we discuss in 
detail the scaling limit construction of the field theory limit given in the previous lecture 
and present the relation of the two spin correlation to Painleve functions. For the inho- 
mogeneous case we discuss the breakdown of this scaling picture which is caused by the 
spreading out of the critical point into a critical region. 

1. Introduction 

The oldest and most completely studied system in statistical mechanics which has 
finite range forces and has a critical point is the two dimensional Ising model in zero 
magnetic field. Consequently, by means of the connection between statistical mechanics 
presented in the previous lecture, the Ising model is also the most extensively studied 
quantum field theory. The principal purpose of this lecture is the presentation of the theory 
of this model. This will concretely illustrate many of the general principles previously 
presented. 

However the most surprising feature of the Ising model is that there are many con- 
nections with quantum field theory that go beyond the theory of the first lecture. The 
Ising model has led to developments in mathematics which have been widely applied to the 
theory of random matrices, to quantum gravity, and the Ising correlations themselves are 
directly related to N = 2 supersymmetric quantum field theory in two dimensions. There- 
fore a second purpose of this lecture will be a to briefly indicate some of these further 
connections between SM and QFT. 

The two dimensional Ising models we will consider in this lecture will be various special 
cases of the general nearest neighbor inhomogeneous model in a magnetic field specified 
by the interaction energy 



where a^k = ±1, j(k) specifies the row (column) of the lattice, periodic boundary condi- 
tions are assumed and in the most general inhomogeneous case the interaction constants 
E v (j,k) and E h (j,k) are allowed to vary from site to site. We will be particularly con- 
cerned with 1) the homogeneous case where E v and E h are independent of position and 




(1.1) 



i=i k=l 
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2) the case of layered impurities where E v and E h depend only on the row index m and 
H = 0. 

For the case of the homogeneous Ising model the discussion of the previous lecture 
has given an unambiguous prescription to obtain a field theory by constructing a scaling 
limit. The primary goal of this lecture is to carry out this construction in detail for the 
Ising model. In sec. 2 we summarize briefly the various methods used to solve the Ising 
model at H = on the lattice and in sec. 3 we present the results of these studies. In 
sec. 4 we use these results to construct the field theory limit and to present the important 
relations the spin correlations have with Painleve functions. We also indicate how these 
considerations arise in other field theory contexts and we discuss the particle content of 
the theory. We conclude in sec. 5 with a discussion of the inhomogeneous model and the 
new physics which emerges in this situation. 

The Ising model was first solved in one dimension by E. Ising in 1925 and for that 
reason it now bears his name. Since then its study has involved some of the most creative 
developments in mathematical physics. The papers on this model are many starting with 
the original work of Ising and continuing to the present day. I have attempted to list 
the papers most relevant to the present considerations as items pj|- ||137|| of the reference 
section. 

In any bibliography there is inevitably some selection criteria that is applied and I 
should make explicit the criteria used here. Firstly, with possibly the exception of []56|| , I 
have included only original material and excluded the large amount of pedagogic writing 
on the subject. Secondly I have excluded the many applications of Ising model theory to 
physical systems. 

There is one further selection that has been made which is more subtle and needs a bit 
of explanation. There are two important relations between the homogeneous Ising model 
and the XY spin chain in a magnetic field defined by 

Hxy = - ^2 (1 + 7)«+i + \iX ~ 7)«+i + ha z m ) (1.2) 

m 

where cr x , a y and a z are the three Pauli spin matrices. Both of these relations involve a 
commutation relation of Hxy with a transfer matrix of the Ising model. 

The first relation is obtained if we consider a row to row transfer matrix of sec. 7 of 
part I which may be explicitly written in terms of 

T- = (2sinh2^)^/ 2 exp(^*^^) and T h = exp(K h £ a^+x) (1-3) 

k=l k=l 

where K v > h = E v > h /kT and tanhif"* = e~ 2KV as 

T R " = (r v ) 1 / 2 T h (T v ) 1 /' 2 . (1.4) 
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Then [0 

[T R '\H XY }=0 (1.5) 

if 1 

cosh2K u * = - and tanh2K /l = (1 - 7) 1/2 //i. (1-6) 

7 

From this commutation relation it follows that the correlation of two Ising spins in the 
same row are given in terms of XY correlations as 

< o n . m o n , m , > = cosh 2 K°* < a x m a x m , > X y - sinh 2 K v * < a y m o v m , > X y ■ (1-7) 

The second relation is obtained if we consider the diagonal transfer matrix defined by 
(1.7.15) as 

T { % W} = Ue^e Kh ^ + K (1.8) 

fc=i 

This commutes with the special case of the XY Hamiltonian ( |1.2| ) with 

7 = 1 and h = (sinh21f /l sinh 2K V )~ 1 (1.9) 

and thus we find that the correlation of two spins on a diagonal in the Ising model are 
given in terms of the XY correlations as 

< ^0,0^n,n > = < « >XY ■ (1-10) 

Because of this intimate relation of the XY spin chain with the two-dimensional Ising 
model there are some results for the XY model which will also give results for the Ising 
model. In this lecture we make no attempt to be complete in referencing the literature of 
the XY model and give only those papers for which there is no analogous Ising result in the 
literature. We will return to the XY model particularly considered as a dynamical system 
at finite temperature in the fourth lecture and will give further XY references there. 



2. Methods of Solution for H = 



In the 50 years since Onsager |4]] first computed the free energy, the Ising model at 
H = has been solved by many different methods. Some methods are powerful enough to 
compute correlation functions and some can be generalized to solve many other models in 
two dimensions but there is no method known at the present time which is both general 
and powerful enough to compute all correlation functions of all models. Because our focus 
here is on the Ising model itself we will naturally discuss the most powerful of the methods. 
However for the purpose of orientation we will begin our considerations with a brief sketch 
of five of the principal methods of solution in their historical order: 

1) Onsager's algebra || 

2) Fermionization || 0] 

3) Combinatorial [15j Jl8| 



19 2C 21 



4) Star triangle equation and functional equations 

5) The 399 Solution §7[ 



112 
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2.1. Onsager's algebra 



Onsager in his original 1944 paper begins his considerations with the XY model 

H X y =A + hA x 



Hamiltonian (|1.9| ) written as 



where 



A 



and A x = - °"r 



2.1 



(2.2) 



m=l m=l 

and demonstrates that these two operators are part of a much larger algebra 

[At, Am] =4Gi- m 

[Gi, A m ] =2A m+ i — 2A m _i 

[Gi, G m \ =0 

He also finds the relation which incorporates the fact that the lattice has sites 



(2.3) 



G 



m+Lh 



-CG T 



-G m C 



(2.4) 



where C is the spin reversal operator and from these relations alone is able to compute 
all the eigenvalues of Hxy and the spectrum of the transfer matrix T v T h which is seen 
to be of the form e a oA 0e a 1 A 1 ^ This algebra and method of solution remained unused and 
unexplored until it was found ||138|| in 1985 to be related to the special case of the chiral 
Potts model which is now called superintegrable which is a generalization of ( |2.2[ ) from 
two to N states per site. The algebra follows from the two special cases ||139|| [[140 1 



[A Ql [A Ql [A 0l A 1 ]]] = 16[A 0l A 1 ] and [A u [ y A x , [A^Aq]]] = 1G i A i, M (2.5) 



In greatest generality it has been shown ||141|| | |142|| that the algebra ( |2.3|) is sufficient 
to guarantee that all the eigenvalues of the Hamiltonian are in sets containing (2s) : > max 
elements of the form 

E = a + hb + V^m^l + h 2 + haj) 1/2 where mj = -s, -s + 1, • • • , +s (2.6) 



where s is either integer or half integer and j max and the j max numbers aj vary from set 
to set. So far the only known representations of the algebra have s = 1/2. 
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2.2. Fermionization 

In order to compute correlations in the Ising model the algebra of Onsager does not 
seem to be enough. To compute these correlations we need the exceptionally powerful 
method of fermionization introduced by Kaufman || in 1949. This method is so powerful 
that, unlike Onsager's original method, it can be used to obtain results on the completely 
inhomogeneous lattice flOD with if = and E v (j,k) and E h (j,k) allowed to vary in an 
arbitrary random fashion from point to point. The Ising model is the only model which 
has been studied in this enormous degree of generality. 

The transfer matrix ( |1.3|) written in terms of o~\ is not in a fermionic form because 
the o~\ considered as operators commute on different sites whereas fermions anticommute. 
The key to the fermionization of the Ising model by Kaufman || is the introduction of 
anticommuting operators by means of the Jordan- Wigner transformation 

fc-i fc-i 
^ = 2- 1 / 2 (n^)72i-i, al = 2~ 1 ' 2 (\[P l ) l2 ^ al = -i l2k . ll2k = --P k (2.7) 
l=i i=i 

where 

7fc7z + 7*7fc = (2.8) 
If we write the completely inhomogeneous transfer matrices 

T^ = eMj^K h (j,k)atat +1 ) and TJ +1/2 = exp(^ K v *(j, k)a z k ) (2.9) 
k=l fc=i 

in terms of the 7^, the resulting expressions are still quadratic in the 7^. Consequently 
the diagonalizing of the transfer matrix is reduced to the study of free fermions on an 
(inhomogeneous ) lattice. Furthermore since the local correlations of Ising spins are now 
expressible in terms of the correlation of Jordan- Wigner strings of free Fermi operators and 
since all free Fermi correlations can be reduced by Wick's theorem to Pfaffians of two point 
correlations, it is possible in principle to study all of the correlations even on completely 
inhomogeneous lattices. 

2.3. Combinatorial 

The combinatorial method differs from the two preceding methods in that it does not 
make use of a transfer matrix and treats the vertical and horizontal interactions on an 
equal footing. It starts with the observation that because cr 2 = 1 we have the identity 

e Kaa ' = cosh K + aa 1 sinh K. (2.10) 
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which may be used to write the partition function as 

z = ^ e E^=iEi=J^c?.fc)^,k^+i,*+^ h C3,*)^,fc^.*+i] 

= Hl[(coshK v (j,k)co S hK h (j,k)) (2 . n) 

j=ik=i 

II II ^ + fc )c r j,fcCT J + i i fc)(l + ^(j, fc)cr i>fe fTj ifc+ i) 

{cr} j = lk = l 

with ^ = tanhi^j, k). The product on the RHS of (EUlp may be expanded as a sum 
of graphs on the square lattice where on each bond either the factor 1 or ziaa' is used. 
The sum over {a} thus reduces the partition function to the sum of all closed graphs with 
weights of Zi(j : k) for each occupied bond. The combinatorial method then shows that this 
sum may be expressed as the Pfaffian of an antisymmetric matrix and when the lattice is 
homogeneous this Pfaffian is easily computed by Fourier transformation. 

Correlation functions are treated in a similar fashion. We start with the observation 
that again since a 2 = 1 any correlation function such as < a m ,n&m' ,n> > ma y be written 
as 

< ffm,nffm',n' >=< ] [ aa' > (2-12) 

path 

where aa' are nearest neighbors and the path is any path which connects the sites (to, n) 
and (to', n'). Then using ( |2.10| ) we write each factor on the path as 

aa'e Klaa ' = sinh K* + aa' cosh K* = sinhlT(l + z^aa') (2.13) 

and thus correlations are expressed as quotients of a partition function of a lattice with 
the bonds Zi on the path replaced by z~ x divided by the original partition function. The 
numerator partition function is expressed as a Pfaffian exactly as was the original partition 
function and hence all correlation functions can in principle be reduced to the evaluation 
of determinants. 

Even though no explicit use has been made of fermions, operators, or transfer matrices 
this method is actually very close to the fermionization method and if proper choice of the 
arbitrary path is made the two methods can be shown to coincide. 

2.4- Commuting transfer matrices and functional equations 



The diagonal transfer matrix (|1.8|) possesses the very remarkable property that if we 
consider two such transfer matrices T D {u) and T D (u') which are obtained from ( |1.8|) by 
allowing K v and K h to depend on a variable u such that the h of ( |1.9| ) is held fixed that 



[T D (u),T D (u')] = 0. (2.14) 
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This relation seems to have been known to Onsager and follows from a relation on the 
Boltzmann weights [ 



112|| known as the star-triangle equation that 



£ < 

<7 rf =±l 



-&d (Li <J a +L 2 & b +L 3 a c ) 



Re 



-(K 1 a b a c +K 2 cr c cr a +K 3 (T a (T b ) 



if 



sinh 2Kj sinh 2Lj 



k 1 and R 2 = 2k sinh 2L\ sinh 2L2 sinh 2Lj 



(2.15) 



(2.16) 



Using this star triangle equation Baxter [112] demonstrates that the transfer matrix 



satisfies a simple functional equation and he solves that functional equation to obtain all 
the eigenvalues of the transfer matrix for the homogeneous Ising model. This method 
is quite unlike the previous three methods in that it has been applied to a very large 
number of two dimensional models. However it is more restrictive, in a sense, than the 
fermionization and combinatorial method, in that it is not able to deal with the completely 
inhomogeneous lattice. 

This technique has been extended by Baxter to compute the spontaneous magneti- 
zation of the Ising model by means of the introduction of a new object called the corner 
transfer matrix. Moreover, with the exception of some very recent work |143|1 on the XX Z 
chain, these techniques have not yet been extended to the computation of other correla- 
tion functions. The hope is that eventually for the homogeneous lattice all the results of 
the fermionization method can be reproduced by these more general techniques. However 
further discussion of these developments is outside the scope of these lectures. 



2.5. The 399th Solution 



The final method of solution fl87j is an ingenious use of the star triangle equation to 
compute the free energy of the Ising model without recourse to functional equations. This 
method has been extended to compute the free energy of the chiral Potts model [144| but 
thus far it has not been extended to obtain all the other eigenvalues of the transfer matrix. 



3. The homogeneous Ising model at H = 

After this brief discussion of the methods used to obtain results for the Ising model 
at H = we are now able to summarize the results of the computations and use them to 
illustrate the scaling theory and field theory construction outlined in the previous lecture. 
It should be kept in mind, however, that in fact much of the general picture of the first 
lecture was developed after many of these Ising results were found and that the general 
theory was developed as an attempt to understand the Ising model in a much wider context. 
It is fair to say that most of the general theory of critical phenomena originates in the results 
we will present below. 
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3.1. Free energy 

The first analytic result is the famous computation of Onsager [HJ of the free energy 



F = -kT{\n 2 + - (2tt) " 2 / d^i / ln[cosh 2/3£ /l cosh 2/3E" 

— sinh 2f3E h cos #i - sinh 2f3E v cos 2 ] } 



(3.1) 



where we use (3 = 1/kT. This free energy is singular at the unique value of T = T c 
determined from the condition 



sinh 2E v /kT c sinh 2E h /kT c = 1 



(3.2) 



which had been previously found by Kramers and Wannier |[5[ . 

The internal energy u is defined as and from the definition of F is written in 
terms of the nearest neighbor correlations as 



u 



d[3F 



-E h < cto,o0"o,i > -E v < ao,o0i,o > • 



(3.3) 



Hence from (3J.) we obtain 



< 0o,o0o,i >= 



2tt 



2tt 



d.e 



(l-a 1 e lf, )(l-a 2 e 
(1 - aie- i6 )(l - a 2 e i 



with 



cki = e~ 2RV tsmhK h and a 2 = e~ 2RV cothii'' 1 
and < 0o,o0o,i >=< 0o,o0i,o > with K h and K v interchanged. At T c we find 



(3.4) 



(3.5) 



< 00,000,1 >= - coth2(3 c E h gd2(3 c E h 



n 



(3.6) 



where gdx = arctansinhx is the Gudermannian of x. From the internal energy we obtain 
the specific heat and find that near T c 

c= ^ ^ -^^(E h2 smh 2 2p c E v +2E v E h + E v2 smh 2 2/3 c E h )\n\l-T/T c \. (3.7) 
oT it 

Thus the specific heat has a logarithmic divergence at T c which corresponds to the critical 
index a = of the previous lecture. 
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3.2. Transfer matrix spectrum 

Onsager Q also computed the complete excitation spectrum of the transfer matrix. 
He found that the full set of 2 L eigenvalues is composed of 2 sets given as (for L h even) 

r h 

L h 1 
In A+ = — ln(2 sinh 2K h ) + - ^ ±j 2r -i (3.8) 



2 

r=l 



and 

h 



In A" = %- ln(2 sinh 2K h ) ± (K v * - K h ) + \ ±72r (3-9) 

r=l 

where 7 r is the positive solution of 

cosh 7r = cosh 2K V * cosh 2K h - sinh 2K V * sinh 2K h cos ^- . (3.10) 



rir 



The eigenvectors corresponding to A + (A - ) are even (odd) under spin reflection and the 
± are chosen independently with the restriction that the number of — signs is even (odd). 
When T < T c [K h > K v *) the maximum are asymptotically degenerate as L h — > oo 
whereas for T > T c X + is the maximum eigenvalue and is non degenerate. 

3.3. Spontaneous magnetization 

The spontaneous magnetization was announced by Onsager H in 1949 and a derivation 



was given by Yang [|1C§ in 1952. The remarkably simple result is 

M_(T) = (1 - (sinh2if /l sinh2iT ; )- 2 ) 1 / 8 . (3.11) 
Consequently the critical exponent (3 = 1/8. 

3.4- Spin-spin correlations 

The simplest result for the spin-spin correlation is for the diagonal correlation at 
T = T c where we have the remarkably simple expression 

2 N ~ 1 1 

< 0-0,00*,* >= (-) N n [1 - (3.12) 

i=i 

which as iV — *• oo behaves as 

< o- 0)0 a N , N >~ AN' 1 / 4 (3.13) 
where A is a transcendental constant ~ .6450 • ■ ■ . Thus the critical exponent n = 1/4. 



44 



More generally all spin correlation functions may be expressed as finite determi- 
nants |22] |5(J with a size that depends on the separation of the spins. The simplest 



of these determinental expressions are 



with 



and 



with 



< CO,OOo,AT > = 



Oo 

a i 



a-i 



ajsr-i cin-2 



a-N+i 
a-N+2 

a 



1 

2^ 



d9e 



-in9 



1 -aie iy )(l -a 2 e" 
[l-a 1 e- i0 )(l-a 2 e ie ) 



AM 1/2 



a 
ai 



a_i 
a 



Ojv-1 «7V-2 



a-N+2 
So 



1 

2^ 



— / d6>e 



-in9 



sinh 21^ sinh 21^ - e~ ie 
sinh sinh 2K h - e id 



1/2 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



These finite determinental expressions are particularly effective in studying the correlations 
for small separations. 

The study of the large distance expansions of the two point function for T ^ T c is 
significantly more involved than the T = T c expansion. To effectively study this case it 
is useful to recognize that there are not one but many determinental expressions for the 
correlations and in particular there is an expression as an infinite (Fredholm) determi- 



nant [31] [56]. From this starting point a complete discussion of the long distance behavior 



was made |66| in 1976. We here reproduce the results of that study. 
For T < T c the two point function has the representation 



< 0o,o0Af,JV >= M 2 _ exp(- ^ F m!n) 



n=l 



where and 



-An 



F ( m% = (-l) n [2^(l - z 2 h )} 2n (2n)- 1 (2n) 

e —iMcj>2j-i—iN4>2i gj n 



>! • • • 



2n 

n 



2i-i 



A(0 2 j_l, faj) Sin 2 (02j + 02j+2) 



(3.18) 



(3.19) 
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Here 4 n+i = 4>i 04n+2 = ^2,Wj < and 

A(0 2j _i, <j>2j) = (1 + z 2 h )(l + z 2 v ) - 2z v (l - z 2 h ) cos02j-i - 2^(1 - COS02J. (3.20) 
For T > T c the corresponding result is (in the more symmetric form of [97] [|106|| ) 

oo 

< *o,octm,n >= MlX(M, N) exp (- ^ (3.21) 

n=l 

where 

oo 

X(M, N) = Y^ X( 2n+1) (M, N) (3.22) 



n=0 
71" />7r 



X (2n+1) (M5 N) = (27r) -4n-2 T ^ . . . C 

J —IT J —IT 



2 "+ 1 g-i(M-l)^ 2 j-i-t(iV-l)^ 

II A(o 2/ ' (3-23) 

g— i<t>2j-l _ g— i<t>2j+i 

2izv{l - zg){ 

3=1 



2n 



Ff[2^(l - rj—^r )] 

i=i 

with 4 n+i = 0i, 04n+2 = 02, Im0j < and 

M + = [(sinh2K u sinhX /l )- 2 - l] 1 / 4 (3.24) 

If the exponentials in ( |3. 18|) and (|3.21|) are expanded ||108|| these series are seen to 
be the Kallen-Lehmann form factor (soliton) expansions of Part 1 (7.23) where the index 
2n (2n + 1) represents the number of particles in the state. 

The leading terms of the expansion (of the diagonal correlation) for T > T c as N — > oo 

is 

1 k N 

< a 0>0 a N , N >~ (7rAr)1/2(1 _ 4)1/4 + • ■ • ( 3 - 25 ) 

where 

k y = sinh 2K h sinh 2K V . (3.26) 
The correlation length thus behaves near T c as 

C 1 = lnfc> ~ \T-T e \ (3.27) 

and hence the critical exponent v = 1. 

The leading term of the expansion (of the diagonal correlation) for T < T c as N — ► oo 

is 

2k 2(N+1) 

< a ,oa N>N >~ Ml{l + _ 2 - i)2 + • • •} (3.28) 
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where 

k< = (sinh 2K h sinh 2K V )~ 1 . (3.29) 
The correlation length thus behaves near T c as 

2C 1 = 21nfc< ~ 2\T - T c \ (3.30) 

and hence the critical exponent v' = 1. 



4. Scaling theory and Painleve equations. 

We are now in a position to make contact with the scaling theory of part 1. First 
we note that the two scaling laws 2(3 = v'(d — 2 + rj) and dv' = 2 — a' are satisfied by 
a' = 0, /5 = 1/8, r\ = 1/4 and v' = 1. We also note the predictions from 7' = dv' — 2/3 
that 7' = 7/4 and from 77 = 2 — d(<5 — 1)/(<5 + 1) that 5 = 15. We also note that the critical 
exponents above T c equal those below T c . 

The scaling limit of the full correlation function may now be obtained from the defi- 
nition 

G±(r) = lim M± 2 < a 0fi a M ,N > ■ (4.1) 

scaling 



and the results are |6(| for T < T c 

OO 

G_(r)=exp[-^A 2 "/ 2n (r)] (4.2) 



n=l 



with 



f j^n roc roo 2n ^ 



/2nW = / dy-L 

n 



and for T > T c 



where 



with 



^"Ilt fJiu/a fe + 11^ - !) ( 43 ) 

i=i ^ J i=i 

G+(r) =X(r)G-(r) (4.4) 

00 

X(r) = ^A 2 " +1 x 2n+1 (r) (4.5) 



n=Q 

/OO /"OO 
dyi--- J dy 2n +i 



2n+l 2n 1 n (4.6) 

e yj -n- 1 



n (yj-i)^ n ~y~Ty^~i 



2j 
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and in the above A = 1/tt. 

For large r the leading behavior for T > T c is 

G+(r)~-K (r) (4.7) 

TV 

where K n {r) is the modified Bessel function of order n. The Fourier transform of this is a 
single pole and confirms the existence of a particle in the theory. 
For large r and T < T c the leading behavior is 

G_(r) ~ 1 + 7T — 2 [r 2 \K\ (r) - K 2 (r)] - rK (r)^i(r) + ^ 2 (r)] (4.8) 

The Fourier transform of this is not a single particle pole but rather contains a branch 
cut. This is consistent with the fact that the excitation spectrum of the transfer matrix 
has only even particle excitations for T <T C 

But we are able to prove even more. Not only does this scaled dispersion relation 
hold but it was shown in 1973 []57| [5S| that the scaling function has the remarkable 



property that 

1 f°° 1 

G ± (r) = -[l T r 7 (r/2)]r 7 (r/2)- 1 /2 ex p / dd -^[{l - V 2 ) 2 - (r/) 2 ] (4.9) 

Z Jr/2 4 



where rj(9) satisfies the Painleve III equation [p.45|| ||146 



v'^-ivT-^ + v'-v- 1 (4-10) 

with the boundary conditions 

r){6) ~ 1 - 2XK (26) as 6 -f oo where A = 1/tt. (4.11) 

If we set 

v (e) = with t = 29 (4.12) 

the Painleve III equation reduces to the radial sinh- Gordon equation 

d?iP 1# 1 . 

^ + t^ = 2 Smh2l/ ' (413) 
and the correlation function becomes 

'sinh^(r)^ ( l f°° j r r dip 
l 

2 

In addition if we set 



<ilnG± , „ , „, 

C = r— —± 4.15 
ar 
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that for both G + and G- we have [ 104 



w\2 



K") 2 = 4«' - C) 2 - 4(CTK' - C) + (C) 



(4.16) 



which has an intimate relation with Painleve V [ 147 |- |P-52[ . 

It is of interest to relate these nonlinear differential equations for T ^ T c to the 
previously studied case of T = T c . To do this ||126|| we explicitly introduce a mass scale 
m into the equation ( |4.1ti| ) as r = mf and then let m — > 0. In this limit the nonlinearity 
of ( |4.16[ ) cancels and one integral can be done to yield 



,d 2 G 



+ 2(c + 



l s dG 



dr 2 v " ' A'' dr ' LV ' A' 1 



+ [(c + 



- \}G = 0. 



(4.17) 



If the integration constant c is chosen to be zero and if we write f = (zz*) 1 /2 then, 
considered as an equation in z, the equation ( |4.17| ) is recognized as the second order 
equation satisfied by a second level degenerate operator in conformal field theory ||153|| [ |154 | . 

Moreover, once these nonlinear differential equations have been found for the scaling 
limit it is to be expected that similar results exist for the full correlations on the lattice as 
given by the dispersion relations ( |3.18| ) and ( |3.21| ). In fact there are several such results 
known. 

The first of these generalizations are quadratic difference equations for the correlation 
functions on the lattice. These are most compactly displayed if we note that the correlation 
functions above and below T c are related by what is called a duality relation. Defining the 
new constants K v * and K h * by 

sinh 2K V * sinh 2K h = 1 and sinh 2K h * sinh 2K V = 1 (4.18) 

and letting < 0"o,oCM,jv >* to be the spin-spin correlation on this dual lattice we have 



< cto,oCm,jv > 



(sinh 2K V sinh 2K h ) 1 / 2 X(M, N) . 



Thus we find the two nonlinear partial difference equations [97] |)8[] [|106 | 

2 

< O~0,0O~M.N > — < 0"O,O0"M-1,AT >< CTo,oO"M+l,Ar > 

= - sum" 2 2^ < a 0fi a M ,N >* 2 

— < C"0,0CM,;v-1 >*< cr 0) oCM,AT+l >" 



(4.19) 



(4.20) 



and 



< O"0,0CM,7V > — < Crofi&M,N-l >< O"0,0C"Af,JV+l > 

= - sinh -2 2K h ( < a , a M ,N >* 2 

_ < C0,0O~M-l,N >*< O"0,0OM+l,AT >*)■ 



(4.21) 
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At T = T c we have < cr 0] oO"M,Ar >=< cro,o cr M,iV >* and these two equations reduce to the 
remarkable simple result which is valid except at M = N = 



Sinh21f (< (To$(TM,N > 2 - < (To,O^M-l,N >< (TQ,Q(TM+1,N >) 

+ sinh2i^ /l (< a , a M ,N > 2 - < cro.ocr^jv-i >< Oo,oOm,7V+i >) = 0. 



(4.22) 



This is the discrete (imaginary) time Toda equation ||155|| introduced by Hirota [ |156|| . 
We may also express the lattice correlations solely in terms of X(M, N) which generalizes 
( |4.9| )and derive a partial difference equation for X(M, N) which is a lattice generalization 
of the Painleve equation ( [4. 1U| ) . For details see |lU6| j. 

The second generalization is to nonlinear differential equations in the temperature. 
The simplest example is for the diagonal correlation function < 00 ,o&n,n >• Here we 
define for T < T c 

t = (sinh 2^ sinh 2if 71 ) 2 (4.23) 



and 



and for T > T c 



and 



a(t) = t[t - 1) — In < a 0t0 a NjN > --. 



t = (sinh 2K V sinh 2K } 



/ \ / i\ d , 1 
a(t) = t[t - 1)— In < cr ,ocrAf,Af > —jt 



(4.24) 
(4.25) 
(4.26) 



where in both cases t > 1. Then it was found by Jimbo and Miwa [|104| that 



„. d 2 a. 



, T o r , ^^da l9 da ^^da 



1,, da . 

'-iK'*-')- 



(4.27) 



This is equivalent to the sixth Painleve equation [ [f50|| . For the correlations off the diagonal 
there also exist systems of such differential equations. We refer the reader to ||150|| for 
details. 

But we have not yet succeeded in verifying that the functions G±(r) satisfy all the 
properties of a scaling function. To do this we must show that as r ^ we regain the 
T = T c expression of (|3.13|) and show that 77 = 1/4. This was first done in 1977 when it 
was shown that [ffOl that as r — » 



1 



7](r/2) = Bt a {l - i~;B- 2 (l - arr 
lb 



2„,2-cr 



+ o(t 2 )} 



(4.28) 
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where the two constants B and a are determined in terms of A as 



<r(A) = — arcsin(7rA), B = B(a) 

7T 



_ 3(T r((l-a)/2) 



r((l + a)/2) 



(4.29) 



When A — > 1/tt then a — > 1 and the correlation behaves as 



G±(r) 



const r 1 / 4 {1 ± -r[ln(-r) + 7^] H r 2 + 

2 8 16 



(4.30) 



where 7^ is Euler's constant. 

To complete the verification of scaling theory we must show that the constant in 
(fO0|) is the one obtained from the transcendental constant A of (|3.13|) . This was done 
by Tracy ||134|| ||135|| in 1991. Thus for the two point function all the features in the 
definition of the scaling theory have been verified. It took 47 years to achieve this level of 
understanding. 

Similar studies have been made for all the higher multipoint correlations. Dispersion 
relations on the lattice have been found [71] , sets of holononic partial differential equations 
have been found in the scaling limit f73|| , |85|] , |9l|1 -|l93|| , ||102|| , and quadratic difference equa- 



tions exist on the lattice for all temperatures ||107|| . These multipoint correlations satisfy 
the scaling properties and the Osterwalder-Schrader axioms [ 108 ] 1 115|| . Moreover the S 
matrix has been computed for T > T c ]72| . The scattering is completely elastic and in the 
n particle channel the S matrix is (— l) n ( n-1 )/ 2 . Thus both for T above and below T c this 
Ising field theory constructed from statistical mechanics satisfies precisely the properties 
needed for a relativistic field theory. This is the only interacting field theory for which this 
construction has been explicitly carried out. 

These expressions of the Ising model correlations in terms of Painleve functions are 
some of the most beautiful results in mathematical physics. But what is even more re- 
markable is the fact that these are not isolated results but form part of a very much larger 
mathematical structure which, though it was first glimpsed in the Ising model, has a very 
large range of applicability. 

One extension of these results is to the correlation functions of other specific statistical 
and field theory models: the impenetrable Bose | 157|| -[ 168 | and Fermi | 169j| [ 170 ] gases, 
the XY spin chain |98| ||117|| |P-71| - p8^| and the closely related sine-Gordon field theory 
at the decoupling point |p.83[ ||184j| . 

But the mathematical structure has extensions far beyond these specific models. One 
part of this structure is the theory of isomonodromic deformations of linear differential 
equations and their relation to determinants of singular operators. These determinants 
are called r functions and are mathematical generalizations of the determinants of the 
Ising correlations. One presentation of this theory is ||149||— This treatment not 
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only encompasses the six second order Painleve functions but also establishes the general 
principles that derive all higher order equations with the Painleve property. 

Another part of the mathematical structure is the connection which the theory of 
random matrices 1118511 has with Painleve functions. This connection was first found for 



the Gaussian ensembles by Jimbo, Miwa, Mori, and Sato [157] in 1980 and subsequently 
vastly extended to many other ensembles of random matrices |186|| -[ 2"0~0" |. Moreover the 



random matrix problems have been connected with 2 dimensional quantum gravity ||186 



| 191|| and this forms another major connection between statistical mechanics and quantum 
field theory. 

Finally we note that exactly the Ising correlations given by the PHI equation have 
recently been seen to arise in iV = 2 supersymmetric quantum field theories [|201|1 . However 
in this case values of the boundary condition constant A other than 1/tt lead to physical 
results. Further discussion of these connections to quantum field theory are beyond the 
scope of these lectures but it is surely expected that there will be more applications and 
generalizations of the Ising model results. 

It remains to consider the particle interpretation of the Ising field theory. But now 
it may appear as if we have encountered a paradox. On the one hand we have faithfully 
executed all the steps needed to show that the theory has interactions. However, on the 
other hand the actual computations were done by a formalism that related the Ising model 
to a free Fermi system. From that point of view the reader may ask what right do I have 
to call the theory anything but free. 

The resolution of this question is most interesting. We emphasized in the first lecture 
in discussing the concept of particle that we really have not given a definition of a particle as 
a thing in itself but only in its interaction with those fields which we chose to call observable. 
In the interpretation of the Ising model as a spin model we are certainly conceiving of the 
spins being observable by their interactions with the external magnetic field H. Indeed 
neutron scattering experiments do exactly that and the Fourier transform of the two point 
spin correlation is what it is claimed is measured. In terms of this interpretation and 
measurement the computations discussed above are certainly to be interpreted as showing 
the presence of strongly interacting bosonic excitations. The construction of the Fermi 
operators is very non-local in terms of the spins and therefore will not couple to the same 
external fields that the spins do. Indeed it is hard to conceive of a possible experiment 
that will couple to these fermions. Thus the fermion Greens functions which are indeed 
free and which are indeed used in the intermediate stages of some methods of computing 
the Ising spin correlations are not excited by the possible laboratory experiment. The 
observed particles are the interacting bosons and not the free fermions. 

In fact this duality between bosons and fermions seen in these computations is not 
unique to the Ising model and in fact is very common in one dimensional quantum systems 
and the related two dimensional classical systems. Much as it may seem to violate naive 
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expectations one local field theory can often be transformed into another local field theory 
by means of a non local transformation. This may indeed be an effect that cannot occur 
in higher dimensions (although I do not know of a definitive study) but because the effect 
exists we must be very cautious when giving a particle interpretation to a quantum field 
theory. In the end it must never be forgotten that it is the measuring apparatus which is 
the authority that tells us which interpretation is physically observed. 

Finally we remark about the 4 e perturbation theory from d=4. Such an expansion 
gives results for the critical exponents of the scaling theory. This expansion does not 
converge but is hopefully asymptotic and Borel summable. The results in two dimensions 
reproduce the known results with a remarkably good accuracy [[202| . But for field theory 
applications these critical exponent computations are not really what is desired. What we 
really want are the Greens functions. So far there has been little attempt to regain the 
Painleve property of the Greens functions or the S matrix from the e expansion. 

We close with the observation that in the e expansion it is tacitly assumed that the 
particle spectrum in d=4 is the same as in d=3 or d=2. As long as this is the case there 
is some hope that the method can succeed. But if the genuine particle spectra in d=3 or 
d—2 is not the same then it is extremely doubtful that the e perturbation expansion can 
be powerful enough to give the Greens functions of the theory. We will return to this in 
the next lecture since for H ^ the question of how many particles are in the spectrum is 
a matter of major concern. 



5. The layered Ising model at if = 

The results of the previous sections have been restricted to the homogeneous Ising 
model. However the fermionization and combinatorial methods are not restricted to 
this case and their initial formulation is valid also for the completely inhomogeneous lat- 
tice ( p. . 1|) where the bonds E v (j, k) and E h (j, k) are completely arbitrary. For this reason 
many of the difference equations satisfied by the correlations are valid for the completely 



random lattice [107]. However, in order to get complete results it has been necessary to 



impose translational invariance in one direction. Therefore in this final section we will 
consider the layered Ising model where E v (j,k) = E v (j). We will also consider for con- 
venience E h (j 1 k) = E but this specialization can easily be removed and does not affect 
any of the physics. 

We begin our considerations with the case where the N bonds E v (j) for 1 < j < N 
are arbitrarily chosen but then this configuration is repeated to make up the full lattice. 
Thus in this case the translational invariance in the vertical direction is reduced but not 
totally destroyed. The free energy of this model has been studied ]6(| where it was shown 
that there is still a unique critical temperature determined from 

N 

i=i 
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At this critical temperature the specific heat still has a logarithmic divergence but now 
the amplitude depends on all the iV energies E v (j). Moreover the amplitude depends 
sensitively on the geometric ordering of the N bonds. In general all AH permutations 
of the iV distinct bonds give different amplitudes and furthermore these amplitudes are 
less than what would be obtained is we considered replacing all the vertical bonds by the 
average value. In this way we argue that there is a sense in which the logarithmic divergence 
can be said to be caused by all the bonds acting together in a cooperative fashion. The 
question of greatest interest then whether the destruction of translational invariance can 
also destroy the logarithmic singularity by driving the amplitude of the logarithm to zero. 

There are of course many ways to destroy translational invariance. The mildest way 
is to have the vertical bonds form a sort of quasi crystal. A particularly nice example of 



this is the Fibonacci lattice considered by Tracy [ 129 1 . This lattice is defined as follows 



Consider a set of sequences S n of the letters A and B defined recursively by 

S n+1 = S n S n -i, with S = B, Si = A (5.2) 

For example S2 = AB, S3 = ABA, and S4 = ABAAB. The sequence S n contains F n _i 
A's where F n are the Fibonacci numbers defined by 

F n+1 = F n + F n _!, F = Fx = l (5.3) 

Now consider two values of the energies E v = E\ and E v Bl place them in the lattice 
according to the sequence S n and repeat the configuration to build up the complete lattice. 
In the limit that the sequence length goes to infinity, so that the translational invariance 
is destroyed, it is found that there is still a single transition temperature obtained from 



e~ 2Eh/kT < = z: Ac zf Bc (5-4) 

where 



a 



- 1 = (l + 2 1 / 2 )/2 (5.5) 



and at this temperature the specific heat still has a logarithmic divergence with the am- 
plitude 

^ = i^r^C^c - ^ C )l ^^l{2^ + ^«(^"i c - ^a c ) + ^« 2 (^ c - ^ c )} 2 (5.6) 

where x = z v b c / z v a c - Thus it is definitely not the case that the mere destruction of 
translational invariance is sufficient to destroy the logarithmic singularity in the specific 
heat. 
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However if the lattice is made completely random in the sense that each vertical bond 
is chosen as an independent random variable from a probability distribution P(E V ) the 



situation changes drastically jJ6|, jJ5|, @, |E7|], [5? 

The first qualitative change is that now instead of the zeroes of the partition function 
pinching the real temperature axis at a point they pinch in an entire line segment that 
runs from T u to T l where T u ' 1 is the critical temperature that would result if all bonds E v 
were replaced by the maximum (minimum) value allowed by the probability distribution 



P(E V ). These two temperatures are called the Griffiths temperatures of the system ||39|| . 
There is still a singularity in the free energy at the temperature we can still denote as T c 
defined by ( |5.1| ) which is now written as 

0= dE v P(E v )ln(zl c e 4Eh / kT °) (5.7) 
Jo 

but now the singularity is not a logarithmic divergence. There is a new temperature scale 
introduced into the problem when T is between the Griffiths temperatures such that near 
this T c we define 

2f °°dEvp(Ev)Hz 2 v e 4Eh/kT ) 

J^dE-P(E-)[lii(z^e 4Eh / kT )] 2 ' [ ' ' 

and find that the singular part of the specific heat is now proportional to 

d^—lnKsi^-il + cP)- 1 ]. (5.9) 

This function does not diverge at 5 = but instead has an infinitely differentiable essential 
singularity. 

There are many other new features that occur in the physics of this random layered 
Ising model. For example the magnetic susceptibility does not diverge at T c but is infinite 
in an entire range of temperatures around T c . The conclusion is that the scaling theory 
developed for the standard second order phase transition with one length scale has totally 
broken down in this random system and that our intuition must be enlarged. 

One interpretation of these computations is that randomness smoothes out singular- 
ities in the specific heat. Therefore, since the one dimensional randomness of the layered 
Ising model has already smoothed out the logarithmic divergence to an infinitely differen- 
tiable essential singularity it can be argued that making the system random in the second 
direction cannot make the specific heat more singular. This would argue that the specific 
heat of the fully random Ising model should be finite. 



This argument has been opposed by Dotsenko and Dotsenko [|203|1 who study a system 
with E v (j, k) and E h (j, k) all chosen as independent random variables with the probability 
distribution 

P(E) = (l-p)8(E-E 1 ) + P 5(E-E 2 ). (5.10) 
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If p is considered small and if only the lowest terms in p are kept then it is found [203] that 
the specific heat has a divergence of the form ln(ln|T — T c | _1 ). However the fact that 
only lowest orders in p are kept means that effectively the distinction between the upper 
and lower Griffiths temperatures is ignored and the transition seems to be replaced by 
some effective transition in which the zeroes still pinch at a point. This would seem to 
exclude from consideration the temperature scale inside the Griffiths temperatures where 
the actual transition takes place. Hence it may be argued that this computation ||203]| does 
not incorporate the new scales seen in the layered Ising model and is not applicable directly 
at the point of singularity itself. 

More recently the fully random Ising model has been considered by Ziegler 



in a field theoretic fashion which does incorporate the new length scales and explicitly 
considers the new behavior between the Griffiths temperatures. These considerations do 



lead to a finite specific heat. Furthermore D. Fisher [136] has applied renormalization 



group arguments to the XY random spin chain related to the layered Ising model and 
finds that the exponent of the spontaneous mag netization is (3 - 5 1 / 2 )/2. There are still 
many open questions in the theory of random Ising models and their study will continue 
to enlarge our intuition of the effects of impurities on phase transitions. 
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Part III. Ising Field Theory at H ^ 



Abstract 

We discuss the relation of the Ising-gauge Ising Higgs theory to the Ising model in a 
magnetic field. We then present the methods used to study the Ising model in a magnetic 
field and summarize all of the known results for the field theory limit. These results are 
then used to illustrate the phenomena of confinement and to discuss the notion of the size 
of a particle. 

1. Introduction 

The solution of the Ising model in a magnetic field has been considered an outstanding 
problem ever since the solution of the zero field free energy 50 years ago. At first it was 
considered because of the relevance to magnetism but in the early 70's it achieved even 
more importance when it was recognized [Q [[| as the simplest non trivial realization of 
the lattice gauge theory of Wilson . 

In the 50 years in which this problem has been investigated there have been two main 
lines of investigation. The first is to use our detailed knowledge of the correlation functions 
at zero magnetic field as the basis of an expansion for small values of H. This leads || to 
the picture that a one parameter family of scaling field theories can be constructed in the 
limit 

H -> 0, T -> T c with h = const H/\T - T c | 15/8 fixed. (1.1) 

The particle content of these field theories depends on the value of h and for T < T c and 
small values of h the spectrum contains a very large number of particles which become 
dense as h — > 0. 

Much more recently it has been discovered by Zamolodchikov that at T = T c 
(h = oo) the field theory limit has the remarkable property that it is integrable. The 
spectrum has been computed to contain 8 particles and the S matrices of the scattering of 
these particles has been found. 

In this lecture we will present these results. In sec. 2 we will present the connection of 
the Ising model and lattice gauge theory. In sec. 3 we will sketch the method of expansion 
for small H. In sec. 4 we will sketch the method of computations for T = T c . In sec. 5 we 
summarize the known results and we conclude in sec. 6 by using these results to discuss 
the physics of confinement. We conclude in sec. 5 by using these results to discuss the 
physics of confinement. 
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2. Ising— Gauge Ising— Higgs 

In (5.6) of the first lecture we introduced the general lattice gauge theory of Wil- 
son |J for a single scalar field coupled to a gauge field. However in the subsequent 
discussion this action was never really discussed and, in particular, in our discussion of 
particles and scattering we made no mention of any particular features connected with 
gauge theory at all. But in applications to particle physics gauge theories are of extreme 
importance. Therefore if we are to have any genuine understanding of the connection be- 
tween statistical mechanics and quantum field theory we must consider gauge theories in 
detail. If we could not do this we would be in the position of having two types of field 
theory-nongauge theories like the n vector model that are relevant to statistical physics and 
gauge theories which are relevant to particle physics. It is thus extremely fortunate that in 
the early 70's it was discovered by Wegner |]J and Balian, Drouffe, and Itzykson that 
there is a connection between the homogeneous Ising model defined by 

E = —E }^ j ((T j,/tcr j+i t k + o j,kC j,k+l) — H 2j cr j >k (2.1) 

j,k j,k 

and the simplest possible lattice gauge theory of Wilson || : the Ising-Gauge Ising-Higgs 
model in two dimensions specified by the Euclidean action 

L 

S = — Eg Sj + l/2,k s j + l,k+l/2 s j + l/2,k+l s j,k+l/2 

7 , A; = 1 . . 

L (2-2) 

-Eh ^2 (0j,k s j + l/2,k<Pj + l,k + (pj,kSj : k+l/20j,k+l)- 
j,k=l 

Here s^ m = ±1, 4>j t k = ±1, j(k) labels the row (column) of the lattice and the variables s^ m 
are associated with the links of the lattice and are specified by the location of the midpoint 
of the link. In the language of gauge theory <f>j^ is called a Higgs field and Sj ^+1/2 is a 
gauge field. We here will show the relation of ( |2.2| ) to (|2.1| ) following closely ref 0. 

The action ( |2.2[ ) has the local gauge invariance property that if we consider any 
particular lattice site (j, k) and replace the four variables s that lie on the links that end 
on (J, k) and the variable <pj t k by their negatives the action remains unchanged. Therefore 
in the computation of expectation values from the general formula 

< 0j >= Z-i E Oje- s / kT , Z = J2J2 e ~ S/kT ( 2 - 3 ) 

only operators Oj invariant under the gauge transformation give non vanishing expectation 
values. We will be particularly interested in the correlations 

< 0j,kSj,k+l/2 ( j ) j,k+l <j>j' ,k'Sj' ,k' + l/2 ( j ) j' ,k+l/2 > (2-4) 
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which we abbreviate as < Hj^Hji^' > and call a Higgs-Higgs correlation and 



< 8j+i/2,kSj+l,k+l/28j+i/2,k+lSj,k+l/2 

(2.5) 

s j' + l/2,k' s j' + l,k' + l/2 s j' + l/2,k' + l s j' ,fe+l/2 > 

which we abbreviate as < Pj^Pj'k' > an d call a plaquette-plaquette correlation. We 
note that the correlations of two gauge variables < Sj ] e+ i/2 s j',k'+i/2 > an d of two Higgs 
variables < 4>j,k<Pj' ,k' > vanish identically because they are not gauge invariant. 

Consider first the partition function Z of (|2.3|) . Because of the local gauge invariance 
the sum over the {0} is trivially done and thus 

Z = 2 L2 J2 e ~ S ' /kT ( 2 - 6 ) 

where S' is obtained from S by setting all tfij^ = 1- Calling Ei/kT = Ki and using the 
fact that Sj &+1/2 and 4>j,k take on only the values ±1 we now write 

e K a ssss = cos hK g (l + z g ssss) and e KhS = coakK h (l + z h s) (2.7) 

where z g ^ = tanhi^/j. Thus 

Z = (2 cosh K g cosh 2 K h ) L2 ^ + z g ssss)(l + z h s) (2.8) 

{s} j,k 

where the indices on s have been suppressed. We may now expand the products and do 
the sum over all s = ±1 and obtain 

Z = (8 cosh Kg cosh 2 K h ) L2 z t z g ( 2 - 9 ) 

graphs 

where the sum is over all closed (intersecting) polygons on the lattice, A is the interior 
area of the polygons and P is the length of the perimeter of the polygons. 

To see the relation with ( |2.1| ) define the variable er^ to be minus one if the plaquette 
whose center is at j 1 + 1/2, k + 1/2 is not used in a term in (2^9) and equal to one if it is 
used. Then 

A = g X^ 1 ~ a i> k ^ Pv = \ X^ 1 ~ °** <7 i> fc +i) and P h = - - a j>k a j+1)k ) (2.10) 

j,k " j,k " j,k 

with P = P h + P v . Thus 

Z = (8 cosh Kg cosh 2 K h ) L2 ^ ex P{^ ln ^<? X^ 1 ~ 

<T = ±1 ~ J fc 

! (2-11) 
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which is clearly proportional to the homogeneous Ising model (2J.) with E v (j,k) = 
E h (j, k) = E if 

H/kT= ^nzg and E/kT = -lnz h . (2.12) 

The important point to note in this derivation is that it goes through word for word 
if we let z g and Zh vary from site to site on the lattice. This may now be used to calculate 
correlation functions. Consider, for example < HqHr > . By definition 

< H H R >= 2 L2 Z~ l So,i/2S M ,N+i/2e- S ' /kT . (2.13) 

{*} 

From e~ s ' l kT we extract the factor e Kh8 °' 1 / 2 e Kh8M ' N + 1 / 2 and write 
So,i/2S M ,N+i/2e KhS °' 1/2 e KhSM ' N +^ 

= COsh 2 K h S 0yl/2 S M ,N+l/2(l- + ZhS ,l/2)(l + Z h S M ,N+l/%) ( 2 " 14 ) 
= cosh 2 K h z\{\ + ^ 1 s ,i/ 2 )(l + Z^ 1 S M ,N+l/2) 

where s 2 = 1 has been used. Thus we find 

< H H R >= z 2 g Z'/Z (2.15) 

where the partition function Z' is identical with Z except that the bonds on (0, 1/2) and 
(M,N + 1/2) are replaced by z^ . The argument leading to (|2.11|) may now be made 



in an identical fashion and we find that Z' is given by ( |2.11| ) with the difference that the 
coefficient of (1 — (j-i^crco) an d (1 — &m-i,n&m,n) is \ bizJ 1 instead of \ \nz g . Thus we 
can write 

< H H r >=z\ < exp[-ln^(l - a_i )0 oo,o) exp[-ln^(l - a m -i )n a mtn )} > 

= < exp[ln z h a-i,oa fi] exp[lnz h a M -icrM,N} > 

= < (coshhiZg + sinhln^a_i ) ocr 0i o)(coshln2;/ l + smhln Zh&M-i,N&M,N) > ■ 

(2.16) 

Therefore we find the desired result 

< H H R > - < H > 2 = smh 2 (2E/kT)(< cr- lj0 a ,o^M-i,N^M,N > - < cr- lj0 a ,o > 2 )- 

(2.17) 

where the correlations on the right are for the Ising model specified by ( j2.12|) . 
In an identical fashion we find that 

< P Pr >-<Po > 2 = sinh 2 (2tf/£;T)(< a 0fi a M ,N >-< a ,o > 2 )- (2.18) 
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The correlation of the plaquettes ( |2.18|) is thus the same as the spin-spin correlations 
considered in general in the previous lecture. But the Higgs-Higgs correlation (|2.17| ) in- 
volves an operator we did not explicitly consider in part 1. This operator aj^j^+i is 
denoted by e^ fc and 

(Tj, k (Tj,k+i = e h j,k ( and °j,k<Tj+i,k = e v j,k) ( 2 - 19 ) 



and is distinguished from Oj^ in that it is even under spin reversal. For the complete 
description of the gauge theory must be considered in some sense to be on an equal 
footing with the original operator <jj^ 

There are many further questions that need to be addressed about the physical inter- 
pretation of these gauge invariant correlations, but we will defer further remarks until sec. 5 
after we have summarized our current understanding of the Ising correlation functions. 



3. Methods of computation for small magnetic field 

The most naive procedure to study the Ising model for H ^ is to use the results 
for the correlation functions at H = discussed in the previous lecture and to expand the 
Greens functions as a series in H. Thus if we denote the interaction energy ( |2.1|) as 

E = E Q -Hj2°j,k (3-1) 
we can study the two point function 

„, (—E /kT+H/kT) y~] . <7j lfe ) s 

< vq,qvm,n >h= Z(H) 2_^ao,o(JM,Ne (3.2) 

a 

in a power series in H as 

< ^0,0^Af,JV >h —M{H) 2 

= ^2^}( H / kT ) n ^2 < a G,0 (7 M,NCfM 1 ,N 1 ■■■(TM n ,N n >H=0 



(3.3) 



n=0 Mi,Ni 



where the superscript c indicates that only the connected part of the correlation is used. 
To be specific consider T < T c . Then the correlation functions at H = have the repre- 
sentation 

< ^o,o^M,N(r mi , Nl ■ ■ ■ <J Mn , Nn >= M* +n e Fn + 2 (3.4) 

where M_ is the spontaneous magnetization. We may now take the scaling limit by using 
for each of the coordinates in ( |3.3| ) the scaling variables of part 1 

m = A(|)- 1 M|T - T c \ and n = ^(O^iVlT - T c \. (3.5) 
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Then if 

HM_(T)A(-)A(0)/\T -T c \ 2 ~ const H\T-T C \~ 15/8 = h (3.6) 
is held fixed the scaled Greens function is given as 

oo 

G|(0, r; A) = ~ hn / rfn - • ■ ^G^ +2 (0, f, r\ • • ■ r n ). (3.7) 

n=0 n ' 



The scaling parameter h of (|3.6|) is the parameter referred to in (|1 . 1|) . It can vary from 



to oo and plays the role of an adjustable coupling constant in the gauge theory. 

The stable particles of the theory are obtained as the poles in the two point function 
which are obtained in coordinate space for large r as 

G c 2 (r) = ^2a l (h)K [(2 + k l (h))r)] ~ Tr 1 /^- 1 ^^ Gz(/l)e -^ (3 . 8) 
I i 

We thus can compute the particle spectrum for small h by using the large r expansion of 
the multi point Greens functions in the cluster expansion. Orient the vector r to point 
along the y axis. Then the multi point Greens functions have a remarkable string property 
that the correlation is concentrated about the y axis in the following form 

G£ +a (0, n, r- 2 e" 2 ^ n+2 (^, ^) (3.9) 

where rrti/r and n^/r 1 / 2 are of order one. Inserting this form into ( |3.7|) we see that each 
integral over fi gives a factor of r 3//2 and hence for large r 

GUr; h) ~ r- 2 e~ 2r V ^(r 3 / 2 h) n c n (3.10) 



n U 



where 



o 



c n = {drrii} / {dni}h n+ 2(mi,rii). (3.11) 



Thus we see that r always appears in the combination /ir 3 / 2 . Therefore the only way for 
the form ( |3.8|) to be possible is if 

ki(h)=h 2/3 Xi and a t (h) = ha h (3.12) 

To determine the numbers A; and a/ the detailed form for the scaling functions must be 
used. When this is done it is found that ai is independent of / and that A; is obtained as 
the solutions of 

A 3/2 A 3/2 
</l/3(V) + J -V3(V) =0 ( 3 - 13 ) 
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where J n {z) is the Bessel function of order n. 

This expansion method can also be applied to correlations for spin for T > T c and 
again the scaling parameter h of ( |1.1| ) is a one parameter adjustable coupling constant 
for the gauge theory. Thus to be complete in addition to the value of h we need to specify 
whether the scaling theory is constructed from above or below T c . From this discussion 
might appear that there could be a singularity at h = oo (T = T c , H > 0) since it can be 
approached from two different directions. However the theorem of Lee and Yang || says 
that the only singularities in the real H — T plane are at H = for T < T c . Consequently 
the field theory limits must join analytically at h = oo 

Finally the expansion method can be used to study the correlations for e_j- £. For this 
case, of course, the correlations at H = of arbitrary numbers of cr's and e's must be first 
evaluated by the methods of the second lecture. We will defer the presentation of results 
to section 5 where we summarize all results needed to discuss confinement. 



4. Method of computation for T = T c and H ^ 

When T = T c and H = the methods of the preceding section break down completely 
because h = oo. It is thus extremely useful that for this case it was found in 1988 by 
Zamolodchikov|| |j) that results may be obtained by a completely different method which 
is quite unlike any method I have discussed previously in these lectures. 

In principle you might try to extend the perturbation techniques of the previous 
section to T = T c by using the exact T = T C correlation functions in the cluster expansion. 
Furthermore you might pass to the continuum limit by using the large R form of the 
correlations and symbolically writing 

E Ising = E t =t c -H J a(x)d 2 x. (4.1) 

Unfortunately the integrals in the cluster expansion are badly divergent and without fur- 
ther definition the procedure is meaningless. 

However, the basic reason that the Ising model is solvable at H = is that it has an 
infinite number of conservation laws which we exploited when we related the Ising model 
to free fermions. What Zamolodchikov recognized |] is that if T = T c and if the 
continuum limit is taken then there is a sense in which (|4.1|) preserves some, but not all, 
of the conservation laws of H = for H ^ 0. 

This observation is both remarkable and mysterious. The mystery is that no one has 
found these conservation laws for the original lattice definition ( |2.1| ) of the Ising model 
in the sense that no one has found a family of commuting transfer matrices T(u) which 
contains the Ising model with H ^ 0. Apparently Zamolodchikov has discovered that there 
is some more general concept of partial integrability that does not demand that all states 
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on the finite lattice be of the desired form. He has managed to throw away those states 
which do not obey the conservation laws. 

The second departure Zamolodchikov makes from our previous methods is that he 
implements the conservation laws by directly assuming that the states of the system have 
a particle interpretation in Minkowski space and that these particle states scatter with an 
S matrix that depends on the rapidity variable 8k defined by 

p° + p 1 = m k e dk , p° - p 1 = m k e~ dk (4.2) 

where p° is the energy, p 1 is the momentum, and nik is the mass of the kth particle. He 
then supplements the basic equations of crossing 

S ab (8) = S ab (i7r-8) (4.3) 

and unitarity 

S ab (8)S ab (-8) = 1 (4.4) 

with the equations for a factorizible S matrix and the requirement that the poles in the S 
matrix elements be at the positions determined by the masses of the other particles in the 
problem. Then when the S matrix elements are constrained by the conservation laws the 
S matrix can be computed and it is found that there are 8 particles with the masses 

TV TV 

m± — m, 7«2 = 2mcos— , m^ = 2m cos — , 

5 30 

777,4 =2m2 cos — , m5 = 2m2 cos — , ma = 2m2 cos — , (4.5) 
30' 15 30 V ; 

TV 7tv TV 2tv 

m-7 =4m2 cos — cos — , ma = 4m,2 cos — cos — . 

5 30 5 15 

These masses have the property that they are given in terms of the components Si of the 
Perron- Frobenius vector of the Cartan matrix of the Lie algebra E$ as rrii/mj = Si/Sj 

More recently a lattice model has been found that also seems to be in the same 
universality class as the critical Ising model in a magnetic field. It is called the dilute A3 
model and the variable has 3 states per site. The model depends on a parameter we may call 
H and the free energy behaves as H ~ _ff 1 + 1 / 15 as H — > which is precisely the behavior 
of the critical Ising model. Furthermore the thermodynamics of the associated quantum 
spin chain have been computed very recently |TIJ and from this the mass spectrum ( |4.5| ). 

From these computations it may be said that the mass spectrum of the Ising field 
theory at T = T c has been computed. Clearly it is most desirable to obtain this result 
directly for the original Ising interaction instead of using this indirect argument. 

Finally it is most interesting to remark that only the masses mi, mi and 777.3 are less 
than 2m2- Therefore it might have been expected that the particles of mass m^ — m% would 
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be unstable. This instability cannot occur in an integrable model denned by commuting 
transfer matrices since it is rigorously known that production and decay are forbidden by 
the infinite conservation laws. However it seems highly likely that any deviation of h from 
oo will destroy the infinite number of conservation laws and hence the particles lying above 
the threshold of decay into 2 particles of mass mi should become unstable. It is thus very 
tempting to see if there is some way to perturb this model in this small non integrable 
fashion to obtain the first model computation of particle decay in a 2 dimensional model 
system. 



5. Summary of results 

We may now summarize the results known for the scaled Ising model in a magnetic 
field. The most basic of these results in given in graphical form in Fig. 1 where we 
schematically indicate the masses and threshold singularities for the spectrum as we go 
from a) H = 0, T > T c to b) H ~ 0, T > T c to c) H > 0, T = T c to d) H ~ 0, T < T c and 
finally to e) H = 0, T < T c . Clearly as we move along the curve from 1) to 5) we pull more 
and more masses out from the two particle threshold and eventually at H = 0, T < T c we 
get an infinite number of poles which become a cut. 

If the spin spin correlation couples to all the particles then this gives the singularity 
structure of the scaled two point function G±. However there is some indication that 



not all particles couple to the two point function (at least at T = T c ). It is surely important 
to analytical study this phenomenon. 

To complete the summary we need the corresponding results for the density density 
correlation < eo,o e M,7V > — < eo,o > 2 . At if = it is known [|12] for both T > T c and 
T <T C that 

D±(r) = const lim \T - T c |" 2 (< e 0fi e M ,N >-< e ,o > 2 ) = KI{t) - K 2 {r). (5.1) 

scaling 

Thus D(r) couples only to 2 particle states and not to any higher multi particle states. As 
limiting cases we have 

D±(r) ~ — r~ 2 e~ 2r as r — > oo (5.2) 

and 

D±(r) ~ r~ 2 as r — ► 0. (5-3) 



Finally for h small it has been shown |T3[ that for large r 



D± (r, h) ~ const G± (r, h) . (5.4) 

However this equality is only valid for the leading singularities and in general there is no 
reason to suppose that e and a couple to the same particles in the spectrum. 
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6. Confinement in the Ising-gauge theory 

We may now combine the results of the previous sections to discuss confinement in 
this most simple of all gauge theories coupled to matter. In previous talks at this summer 
school the test for confinement was whether a large Wilson loop behaves with an area or 
a perimeter law. Unfortunatly this criterion is only useful for either pure gauge theories 
or for those for which the matter fields are treated as infinitely heavy or static. As soon 
as matter fields are dynamic the large Wilson loops all follow a perimeter law and thus do 
not distinguish between confined and unconfined phases. Consequently we are forced to 
discuss confinement in terms of the scaled plaquette-plaquette correlation which we call 

P±{r, h) = lim smh~ 2 (2H /kT)M± 2 (T)(< P P R > - < P > 2 ) = G±(r, h) (6.1) 

scaling 

and the scaled Higgs-Higgs correlation which we call 

H±(r, h) = lim const sinh" 2 (2E/A;T) \T - T c \~ 2 (< H H r > - < H > 2 ) 

scaling 2^ 

= D±(r, h). 

Consider first the Higgs-Higgs correlation H± (r, h) . At some naive level this operator 
creates a pair of Higgs particles connected by a gauge string at one point and destroys 
them at some other point. The correlation function is thus a measure of the propaga- 
tion. Consider first H = 0, T > T c . Then from the results of sec. 5 on D + (r, 0) we see 
that H + (r, h) couples only to two particle states and not to one particle states. This is 
interpreted as saying that two free particles are created each of mass m(= 1). 

Next consider h > with T > T c . Now from sec. 5 H + (r,h) couples to a single 
particle state of mass m. This is clearly interpreted as saying that the two Higgs particles 
are confined by the gauge string to form a single meson of mass m. This change in 
total mass from 2m to m is a very strong confinement. As we continue around the circle 
towards H = 0, T < T c the confining potential becomes weaker and weaker and more 
excited states (mesons) are produced. The number of confined mesons increases until 
H = 0, T < T c where the weak confining potential disappears altogether and a two particle 
state of unconfined Higgs particles each of mass m is regained. 

In terms of the Ising model magnetic interpretation what is happening is for T < T c 
and weak magnetic field H is that the field creates a weak linear confining potential between 
two kink states. Assuming that the two kink states obey Schroedinger's equation we must 
solve that equation with a linear potential. The eigenvalues of this equation are exactly the 
Bessel function spectrum we obtained by the cluster expansion. Of course in the strongly 
confining region such a simple picture is not possible. 
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This picture of confinement in the weak coupling regime suggests another phenomenon. 
For the low lying mesons in the weak coupling regime the wave function is going to be 
nonzero over a large region (which spreads out to oo as h — > 0). This certainly suggests 
that there will be a size associated with these weakly confined mesons which can become 
quite large. Consider a hypothetical experiment which detects the scattering of the single 
meson in the strongly confining regime. The size grows constantly as h — > for T < T c . 
But for a finite size detector the meson will eventually become larger than the detector 
and hence it will be undetectable. This seems a good illustration of the care that must be 
made in interpreting the results of quantum field theory computations. 

Now consider the other gauge invariant correlation P± (r, h) . This has no Higgs fields 
in it and thus seems to couple not to Higgs but only to the gauge field. The question now 
arises if P± (r, h) is coupled to anything other than mesons and if such excitations could 
possibly be identified with excitations of the gauge field (sometimes called glueballs). If 
indeed there are poles in P± (r, h) which are not poles of H± (r, h) it would seem that there 
are excitations of the gauge field separate from the mesons of confined Higgs particles. 

In conclusion I draw the readers attention to the fact that this gauge theory in- 
terpretation does not seem quite to be the same as the particle interpretation used by 
Zamolodchikov in his bootstrap computation. In order to fully understand scattering in 
Zamolodchikov's interpretation one would like to find in terms of the Ising model operator 
a operators which create each of the eight particles separately. These are presumably the 
interpolating fields which are much discussed in quantum field theory and their construc- 
tion should shed light on the vague concept of the size or extension of particles. 

Figure Captions 

Fig.l. The masses and threshold singularities of the Ising model with a scaled magnetic 
field h. 

a) H = 0, T > T c 

b) H ~ 0, T > T c 

c) H > 0, T = T c 

d) H ~ 0, T < T c 

e) H = 0, T < T c . 
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Part IV. Quantum Statistics and the Chiral Potts Model 



Abstract 



We extend the considerations of the previous lectures from classical to quantum sta- 
tistical mechanics. We discuss diffusion in quantum spin chains and the relation this has to 
random matrices. We conclude with a discussion of the chiral Potts model as an example 
of a system where the physics in Euclidean and Minkowski space can be very different. 

1. Introduction 

Thus far in these lectures I have discussed problems in which classical statistical 
mechanics and quantum field theory are very closely related. This very close relation has 
come about because in the systems considered the models have had the feature that it 
is possible to make the analytic continuation from the Euclidean space where classical 
statistical mechanics is naturally interpreted to Minkowski space where quantum field 
theory is most naturally interpreted. However there are many situations in physics where 
these analytic continuations can not be made and where it is more appropriate to directly 
consider quantum statistical mechanics in Minkowski space to begin with. Consequently 
in this final lecture I want to briefly sketch some of these problems in quantum many body 
statistics. In contrast to the previous lectures I will make no attempt to pretend that 
these remarks will be self contained and I will confine myself as far as possible to sketching 
problems and results to the exclusion of explaining methods. 

In sec. 2 I will return to the question of diffusion raised in the first lecture and explain 
in some detail what is known. In sec. 3 I will briefly discuss what is known about the 
relation of diffusion and the theory of random matrices. I will conclude in sec. 4 by 
discussing the chiral Potts model which was introduced and using it to discuss situations 
where Euclidean and Minkowski physics are different. 

2. Quantum diffusion 

In quantum statistics expectation values are computed from 



where if is a quantum Hamiltonian. This formulation is directly in Minkowski space and 
thus is extremely natural to consider time dependent correlations such a 



< O, >= Z" 1 Ti:0 i e~' ff/fcT . Z = Tre" H / fcT 



(2.1) 



< A(t)B(0) >= Z- 1 Tre- iHt Ae iHt Be- H/kT . 



(2.2) 
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It is clear that the study of these correlation functions involves a knowledge of both the 
spectrum of H and the matrix elements of the operators A and B in the basis of the 
eigenstates of H. 

Let us first consider the spectrum of H. In our previous discussion of quantum field 
theory and statistical mechanics we considered the excitations which are of order one above 
the ground state in the thermodynamic limit. More precisely we considered systems where 
in dimension d the limit 

lim E GS /L d exists (2.3) 

L— >oo 



and also where 



lim (E ex - E GS ) exists. (2.4) 

L— >oo 



If in addition we have the property 

N m a 

hmjE ex -E GS ) = J2 E e «( P D ( 2 - 5 ) 

a=l j a , rules 

where N is the number of species of particles, m a is the number of particles in the state, 
and e a (p) are called the single particle energy levels we say that the spectrum is of the 
quasiparticle form. 

It is generally assumed that (but I know of no proof) that most quantum Hamiltonians 
have this quasi-particle form for the order one excitations. This is usually sufficient to 
study the system in the limit that T — > 0. But if we are interested in quantum statistics at 
T > this quasi-particle discussion of the order one low lying excitations is definitely not 
good enough. In fact for any T > these order one excitations above the ground state, 
even though they form a Hilbert space built on these single particle excitations, are not 
sufficient to describe the physics of the system. Indeed these order one excitations form a 
set of measure zero in the total (non separable) space of all states. In quantum statistics 
we are thus forced immediately to deal with the fact that there may be phenomena in 
nature that are not particle like. The aim of this section is to give some insight into the 
study of these phenomena. This insight will be maximized if we get as far away from the 
zero temperature quasi-particle regime as possible. Therefore we will focus on the opposite 
situation and set T = oo. 

To be concrete we will concentrate on the cases where the most exact information is 
known: 

1) The various special cases of the nearest neighbor XYZ model 

L 

H xyz = - >> r ^., + J y S]S? +1 + J Z S Z 3 S Z J+1 + HS*) (2.6) 
i=i 
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where are the three spin S rotation matrices on the site j and in particular 



S\ = \a\ for S=\ (2.7) 

where trj are the three Pauli spin matrices on the site j. When J x = J y this is called the 
XXZ model and we write 

L L 

H xxz = -2jJ2(S]Sj+i + S]Sy +1 + ASffi +1 ) - (2.8) 

3=1 3=1 

and when J x = J y = J z (A = 1) this is the spin S Heisenberg magnet. 
2) The nearest neighbor plus next nearest neighbor Heisenberg chain 

L L 
H n nn = ~J\ Sj • Sj+1 — J2 Sj • Sj + 2- (2-9) 

3=1 3=1 

The first, and perhaps surprising, fact about these quantum spin chains at T = oo is 
that they do not reduce to a classical limit. Moreover the setting of T = oo has not even 
simplified the problem. Indeed for the XYZ chain at T — with S = 1/2 and H = a 
large number of exact results have been computed since Baxter discovered in 1971 that 
the system is integrable. In contrast at T = oo the only exact results for spin correlations 
are at J z = 0. Even the integr ability of the S = ~ XYZ chain has not been exploited to 
any great extent for T > 0. 



At J z = with S = 1/2 on the other hand the Hamiltonian (|2.6| ) is called the XY 
model and the time dependent correlation functions have been studied in great detail. 
At all temperatures the correlations of a x and a y are known to satisfy partial difference- 
differential equations of the Toda type @ M. Many asymptotic results are known for 



large space and time separations PI~ [PTl - These results all stem, ultimately, from the fact 



that the quasi-particle form for the excitations (|2.5|) is not just valid for the low lying 
excitations in the L — > oo limit but is in fact valid for all finite lattices. Thus the spectrum 
is of the free Fermi form form all temperatures. 

From the large number of results available on the XY model we will restrict our 
attention here to T = oo and to the special case of ( |2.8| ) with A = and H = 0. For this 
case it has been found [§J that 

< <(*K(0) >= [J, m _ n |(2Jt)] 2 (2.10) 

where J n (t) is the Bessel function of order n, and ||-|§ 

<a x m (t)a x n (0)>=S n , m e- j2t2 (2.11) 
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The correlation function C 2 (n,t) =< cxQ(0)cr 2 ((t) > satisfies the conservation law 

oo 

C z (n,t) = l (2.12) 



n= — oo 



and therefore can be compared with the spin diffusion forms of part 1. From the exact 
result ( |2.1Q ) we find as t — > oo 

C z (0,t) = cos 2 (2 Jt- j) + ■■■ (2.13) 

and we find that the exact result for the spatial variance cr(t) 2 is 

oo oo 

a 2 (t) = C z (n,t)n 2 = [J n (2Jt)] 2 n 2 = 2(Jt) 2 . (2.14) 

n — oo n— oo 

Both of these results are in contradiction with the spin diffusion forms for t — > oo of 

C z (0,t) ~ (47rDt)- 1 / 2 , a 2 (t)~2Dt. (2.15) 

It is of course not particularly surprising that a system whose spectrum is exactly that 
of free fermions does not have a diffusive character at infinite temperature. The more inter- 
esting question is what is the behavior of C z (n, t) for the spin 1/2 XXZ model ( |2.8| ) (with 
H = 0). Here for A^O the energy spectrum is definitely not that of a free Fermi system 
and the excitations do scatter. On the other hand this system is integrable and the n 
body scattering is factorizible into a series of 2 body scatterings. Therefore any initial 
momentum distribution of n excitations will be preserved in time. This would lead to the 
conjecture that even though the system is not free diffusion still will not occur. 

It is of interest to confirm or deny this conjecture of no spin diffusion in the spin 1/2 
XXZ model on the basis of an exact computation. Indeed for the isotropic Heisenberg case 
A = 1 this problem has been investigated for the last 20 years [|12|-|14|J. However no one 



has found a way to turn the intuition of the argument given above into a computation. 
Instead what has been done is to expand the correlation function into a series in t as 

00 f-iV 
i=o [Al> - 

where by directly expanding the exponentials in (|2.2| ) we have 

Mi l (n)= lim 2- L T r {[[H,[H,---,[H,a z ]]---]]a n } (2.17) 

L^oo v J / 

21 times 
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and we make use of the identity 



Tr{[A, [A,-.-,[A,B}---]]C} 



21 times / 1 _\ 

= (-l)<Tr{([A, [A,--, [A, B] ■ ■ .}})([A, [A, ■ ■ ■ , [A, C] ■ ■ ■]])}. 

V v ' V „ ' 

I times I times 



The longest of these series for the Heisenberg chain [Ti\ have been obtained on the com- 
puter for terms up to t 30 . 

The technical question is now to determine the t — > oo behavior of the infinite system 
from these truncated power series expansions. The most naive thing to do is to just use 
the first 30 terms of ( |2.16| ) as they stand. This, however, does not give a very good ex- 
trapolation and much better results are obtained for the autocorrelation function C z (0, t) 
if we use the exact piece of information that the Fourier transform of this function is non- 
negative. The finite truncation of the series (|2.16| ), however, does not have a nonnegative 
Fourier transform. However there is a lovely theory of Tchebycheff bounds |0J [ 1]J which 



uses this positivity of the Fourier transform and the finite number of coefficients M^(0) 
to produce exact upper and lower bounds for C z (0,t). We give the bounds from [nj in 
Fig. 1. 

It is obvious from these bounds that there are oscillations in the correlation function 
which are reminiscent of the oscillations in the case J z = 0. It is also obvious that the 
asymptotic regime has not yet been reached. Similarly if the spatial variance is stud- 
ied [inj one finds that the asymptotic regime is also not obtained from the first 30 terms. 



Accordingly none of the authors who have studied the problem have been able to make a 
definitive statement about whether or not the spin diffusion form has been confirmed or 
denied. 

These short time series expansions have been extended to several of the other one 
dimensional spin chains which are not integrable in the sense of commuting transfer 



matrices. These systems are the following special cases of ( [2.8| ) with H = and the given 
maximum computed power of t: 

1) S = 1, A = to order t 22 , 

2) S = 1, A = 1 to order t 18 , 

3) S = oo, A = to order t 18 , 

4) S = oo, A = 1 to order t 16 . 

In all of these cases the maximum order is much smaller that the t 30 which was 
computed for the S = 1/2 Heisenberg magnet and as is to be expected these series are 
not long enough to allow any conclusions to be made about their long time behavior. 
Consequently the present status of the short time series studies of the one dimensional spin 
chains is that there is no system for which series are known of sufficient length to determine 
whether or not the spin diffusion forms hold nor can a difference between integrable versus 
non integrable be detected. 
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3. Random Matrices 

We have now discussed all of the results which are available for the time dependent 
correlation functions at infinite temperature for quantum spin models. Thus we have 
exhausted all the results for which a direct comparison with the theory of spin diffusion 
can be made. However, as was stated at the outset these time dependent spin correlation 
functions depend both on the spectrum of H and the matrix elements of the spin operators. 
It is therefore natural to split the problem into two parts and to study the spectrum of 
H by itself and to attempt to find a way to characterize the energy levels which are not 
low lying levels of order one. Indeed it seems most plausible that it should be possible 
to determine whether or not a system has diffusion solely on the basis of the eigenvalue 
spectrum of H. 



Quite recently a very interesting study has been made JI7] of the spin 1/2 Heisenberg 
chain with nearest neighbor and next nearest neighbor interactions ( |2.9| ). For a chain of 
length 20 the eigenvalues have been numerically computed for states of fixed total spin and 
fixed momentum. This leads to a large number of levels for which E max — E m i n = O(L) 
as L — > oo. As noted earlier the levels that have E ex — E m i n and E max — E ex of order 1 



as L — > oo are expected to have the quasi-particle form (|2.5| ) and these are excluded from 
the data. The spacing of the remaining levels is computed and in any region of order one 
the average spacing is normalized to one (where we note that this average spacing varies 
slowly as we go through the entire distribution of eigenvalues.) 

In Fig. 2 we plot the results of such a study for the spin 1/2 chains ([2.8f) and 
( |2.9| ) with 16 sites, momentum P = 27r/16, S z = 1 and H = 0. Here we plot -P(s), the 
probability of the normalized level spacing s which is defined as the true level spacing S 
divided by the local mean level spacing D. We consider the following 4 cases 

a) H nnn for Ji = 1 and J2 = 0.5. This model is not integrable. 

b) H xxz for A = 1. This model is integrable. 

c) H xxz for A = 0.05. This is also integrable. 

d) H xxz for A = 0. This model is not only integrable but the spectrum is that of 
completely free fermions for which the quasi particle property (|2.5| ) holds on the finite 
lattice without the need to take the thermodynamic limit. 

These results have a very remarkable property. For the integrable cases shown in 
Fig. 2b and 2c the normalized level spacing statistics follows the Poisson distribution 

P(s) = e- S . (3.1) 

The probability is maximum at zero spacing and decreases monotonically. But if J2 is 
sufficiently large (J2/J1 = -5) the normalized level spacing statistics of case a) have a 
completely different character. The probability now vanishes at zero level spacing and has 
a single maximum. A more refined characterization is to note that the spacing distribution 
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is well fit by the distribution function of the ensemble of random matrices described by the 
Gaussian Orthogonal Ensemble which is numerically well approximated by the surmise of 



Wigner [18 



TVS , IX 2\ 



p( fl ) = _exp(--0- (3-2) 

(For a further discussion of this random matrix problem the reader is referred to ||19|| ). 
In |17j it is suggested that this difference between Poisson and GOE distinguishes the 
integrable from the non integrable case. However, we see in Fig. 2d for case d) that if 
the spectrum is totally free that the spectrum is not even Poisson but has an enormous 
degeneracy of levels at zero spacing. 



4. Chiral Potts model 

In the discussion of quantum gravity in the first lecture of this series it was emphasized 
that even though there has been much work on Euclidean quantum gravity that in fact 
there are many problems with such a formulation. Furthermore in this summer school 
in the lectures on string theory and quantum gravity the necessity of considering gravity 
directly in Minkowski space has been strongly emphasized. Consequently I wish to conclude 
these lectures on the relation between quantum field theory and statistical mechanics with 
a discussion of a model where the physics in Minkowski space and Euclidean space are 
very different and not connected by a Wick rotation. This is the Chiral Potts model which 
was introduced as a classical statistical mechanical model in Part 1 of this lecture series 
as the classical interaction energy 

L N-l 

e = -J2H {Kiwik+ir + K(°j,k°j+i,k) n } (4.i) 

j,k n=l 

where a^ k = 1. The Boltzmann weights of this model are given in terms of the E v - ,h as 

N-l 

W v > h (n) = exp—Y,E]' h (4.2) 

i=i 

where u = e 27 ™/^ and the transfer matrix is given in terms of these Boltzmann weights as 

i 

T nui'} = n wv ^ - - w ( 4 - 3 ) 

3=1 



The model has been found to be integrable |H|-[[23] if the Boltzmann weights are restricted 



to lie on the manifold [22 
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where a p , b p , c p , d p and a q ,b q ,c q , d q are restricted to lie on the generalized elliptic curve 

a N + Xb N = Xd N , Xa N + b N = X'c N (4.6) 

with 

A' = (1 - A 2 ) 1/2 . (4.7) 
By integrability we mean the commutation relation 

[T p>q ,T M ,} = 0. (4.8) 

The curve ( |4.6| ) has genus g = N 3 — 2N 2 + 1. It is important to note that this integrable 
model will only describe a physical two dimensional classical model if the Boltzmann 
weights are real and positive. This will be the case if the points p and q are restricted to 



a 



■c = u 1/2 bd, |a| = |d|, and \b\ = \c\. (4.9) 



There is also a one dimensional quantum spin chain which corresponds to this two 
dimensional model in a fashion similar to the way the XY model corresponds to the Ising 
model, namely the Hamiltonian introduced by Howes, Kadanoff and den Nijs |20] for the 
special case N = 3 

L N-l 

Hc P = -J2J2 + «n(^j +1 ) n } (4.10) 

j=l n=l 

where we use 

Xj = I N ®--- X •■•®I jV (4.11) 

site j 

Zj = Jjv® ••• Z Jjv, (4.12) 

site j 

In is the N x N identity matrix, the elements of the N x N matrices X and Z are 

X l>m = S hm+1 (mod iV) (4.13) 

and 

Z l>m = 5 hm u l -\ (4.14) 

and 

a n = exp[z(2n - N)<j>/N]/ sm(nn/N) (4.15) 
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a n = Aexp[z(2n - N)<f>/N]/ eaji{im/N). (4.16) 
When the parameters <p : (f> and A are related by 

cos (f> — A cos 4> (4-17) 

the Hamiltonian is obtained from the classical model obtained by considering the limit 
p -> q 

T P:q = 1(1 + const) + uH cp + 0{u 2 ) (4.18) 

where u measures the deviation of p from q. 

The quantum spin Hamiltonian Q4.10|) has several properties: 

1) H cp is hermitian for (p, (f> and A real; 

2) H cp commutes with the spin rotation operator Ylj=i -^-j'l 

3) H cp is translationally invariant; 

4) Even when H cp is hermitian the matrix elements are real only for <fi = <fi = 0. 
From 1) it follows that the eigenvalues of H cp are real and thus that the Hamiltonian 

is a physical energy operator, from 2) it follows that the eigenvalues of the spin rotation 
operator e 2lTl Q/ N for Q = 0, 1, • • • N — 1 are good quantum numbers and from 3) it follows 
that the eigenvalues of the translation operator e lP with P = 2%k / L with k = 0, 1, • • - L — 1 
give the total momentum of the system. These are standard properties which are shared 
with many other of the systems discussed in this set of lectures. However property 4 is 
new and leads to features we have not seen before: 

1) When the matrix elements are complex the system is not time reversal invariant 

2) For <fi 7^ the system is not parity invariant. This is the reason the model is called 
chiral; 

3) The Perron-Frobenius theorem [53] cannot be invoked and thus ground state level 



crossing may occur. 

The Perron-Frobenius theorem has in fact already played a major role in the connec- 
tion between statistical mechanics and quantum field theory. It says that if the off diagonal 
matrix elements of a matrix are positive then the state with the largest eigenvalue is non 
degenerate. For any statistical mechanical system with real interaction energies the Boltz- 
mann weights are clearly real and positive. Therefore the Perron-Frobenius theorem says 
that the eigenstate of the transfer matrix with the maximum eigenvalue (from which we 
obtained the free energy) is never degenerate no matter how many parameters are varied. 

The integrable two dimensional classical model had real positive Boltzmann weights 
when ([Op holds and thus on this manifold the theory of critical phenomena and the rela- 
tion with Euclidean field theory developed in the previous 3 lectures will apply. However 
this manifold is not the same as the manifold where 0, and A are real. Therefore the 
manifold where the Hamiltonian is Hermitian is not the same as the manifold where the 
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statistical system is physical. Thus there is no reason to believe that the physics of the 
Minkowski space quantum Hamiltonian will be related to the Euclidean statistical system. 

In fact the physics of the statistical and quantum system are very different for this 
model. This is most vividly seen in explicit computations done for the special case 
of ( 4.10 ) called superintegrable where <p = <p = n/2. Here it is found that the ground 
state energy is [^5[ |2(J 

N-l , 

EgsIL = -(1 + A) £ F(--, 1; -— — ) (4.19) 



n=l 



where F(a, b; c; z) is the hypergeometric function. This ground state energy is only singular 
at A = 1 and the previous discussions of the connection between statistical mechanics and 
quantum field theory might lead us to expect that the model had a mass gap which vanished 
only at A = 1. 

However this expectation turns out not to be the case. The excitation spectrum of the 
superintegrable case has been explicitly computed. It is found to be of the quasi particle 
form ( |2.5| ) with the fermionic exclusion rule 

P? + P% for j ? k (4.20) 



and for N = 3 there are two distinct types of quasi particle energies [26| [27] which in the 
sector Q = are 



| l-A I 2/3 „ 

3 /' ll - Al , , UJV r u z v r , r 4A 



7T 



e r (n=2\l-X\ + - ' ^(_^ + _^_)[-^--(l-A) 2 ] 1/2 (4-21 



where 



with — oo < v r < cc and 



and 



UjtV r — 1 UJ 2 tV r — 1 t 3 — 1 

e"^ = 1 ~ ^ V (4.22) 

1 — UJV r 

< P r < 2n (4.23) 



e 2s (P )-4[l-A| + -y i dt {v2st)3 - i [__(i_A)] 



(4.24) 



where 

e" lP = V ■ /3 ( 425 ) 
1 - e l7T / 3 v 2s 

with v 2s > and 

— < P 2s < 2tt. (4.26) 
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Furthermore there is a conservation law on the number of excitations m r and m,2 S 



m r + 2m 2s = (mod3). (4.27) 

There are two very important features to be noted in this spectrum. First the mo- 
mentum range QL26| ) is not the full Brillouin zone of 0-2n. This is perhaps not completely 
in accord with the usual intuition about particle excitations but in fact does not violate 
any principle. 

However the most important feature of the excitation energies is the fact the e r (P) is 
not always positive. In fact explicit numerical evaluation of the integral in ( }4.21| ) shows 
that if 

.9013 < A < 1/.9013- •• (4.28) 

then there is a region of momentum P where e r (P) is negative. Thus when ( f4.28| ) holds 
there has been level crossing and the state for which the eigenvalue ( |4.19j ) has been 
computed is no longer the ground state of the system. Therefore a phase transition has 
taken place because the ground state levels have crossed and this transition does not reflect 



itself in a singularity in the ground state energy (|4.19|) . This sort of phase transition is 
outside the scope of the critical phenomena discussed on the first lecture. 

Level crossing transitions such as the one considered here can potentially occur when- 
ever they are not excluded by the Perron-Frobenius theorem. This is precisely the sort 
of physics which the chiral Potts model was originally invented to study ||20|1 . Moreover 
in the chiral Potts model it can be shown that the correlation functions do not have a 
Lorentz invariant form and that space and time occur in a very asymmetric fashion. The 
full study of these effects is beyond the scope of this lecture and, indeed, many of the prop- 
erties which are known for the N = 2 case, which is the Ising model, have not yet been 
computed. However it is the opinion of the author that at least for the super integr able 
case most of the computations done for the Ising model can be extended to this model and 
that when this is done we will have greatly increased our understanding of the physics of 
level crossing transitions. 



5. Conclusion 

I now have come to the end of this lecture series and I hope that I have fulfilled my 
mandate of explaining how statistical mechanics and quantum field theory are connected. 
Indeed, I hope that I have persuaded the reader that far from being the completely sepa- 
rated fields of research which they were considered to be 30 years ago that they are in fact 
different aspects of exactly the same subject. You might try to make a distinction by say- 
ing that statistical mechanics is more general because it allows the explicit consideration 
of vacua which are not rotationally or Lorentz invariant and it allows for diffusion. But on 
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the other hand since cosmology deals with elevated temperatures and compactified string 
theory certainly breaks Lorentz invariance of the underlying 10 or 26 dimensional space 
even this distinction seems to be only of marginal importance. 

It is very impressive that modern day mathematical physics embraces fields as far 
separated as algebraic geometry and the theory of critical phenomena; that gauge theory 
and the theory of magnets are related: that random matrices and two dimensional quantum 
gravity lead to the same type of nonlinear differential equations that characterized the Ising 
model correlation functions. In the beginning of this century some of the deepest insights 
into the emerging science of quantum mechanics came from statistical mechanics. I hope 
that in this series of lectures I have demonstrated that at the end of the century this 
connection of the two subjects continues to provide profound insights into the deepest 
problems of physics. 
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Figure Captions 

Fig. 1. The bounds on the autocorrelation function < o"o(i)crg(0) > at T = oo of the 
Heisenberg spin chain as obtained by Bohm and Leschke in ||14j| . 

Fig. 2 Plot of the histogram of the probability of the normalized level spacing s = S/D 
for various spin chains with 16 sites, momentum P = j| and < a z >= 1. The following 
cases are considered: 

a) H nnn (|2.9|) with J\ = 1 and J2 = 0.5. The dashed line is the Poisson distribu- 
tion (|3.1| ) and the dotted line is the distribution function of the Gaussian Orthogonal 
Ensemble Q. 

b) H xxz ( |2.8| ) with A = 1 and H = 0. The dashed line is the Poisson distribution ( |3.1|) . 

c) H xxz (|2.8| ) with A = 0.05 and H = 0. The dashed line is the Poisson distribution 

(B- 

d) H xxz ( |2.8j ) with A = and H = 0. The dashed line is the Poisson distribution 

(B- 
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